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ABSTRACT 

We  present  infrared  T-band  (3-4  gm)  nuclear  spectra  of  a  large  sample  of  nearby  ultraluminous  infrared  galaxies 
( ULIRGs).  ULIRGs  classified  optically  as  non-Seyfert  galaxies  ( LINERs,  H  n  regions,  and  unclassified)  are  our  main 
targets.  Using  the  3.3  fim  polycyclic  aromatic  hydrocarbon  (PAH)  emission  and  absorption  features  at  3. 1  //m  due  to 
ice-covered  dust  and  at  3.4  jim  produced  by  bare  carbonaceous  dust,  we  search  for  signatures  of  powerful  AGNs 
deeply  buried  along  virtually  all  lines  of  sight.  The  3.3  /im  PAH  emission,  the  signatures  of  starbursts,  is  detected  in 
all  but  two  non-Seyfert  ULIRGs,  but  the  estimated  starburst  magnitudes  can  account  for  only  a  small  fraction  of 
the  infrared  luminosities.  Three  LINER  ULIRGs  show  spectra  typical  of  almost  pure  buried  AGNs,  namely,  strong 
absorption  features  with  very  small  equivalent  width  PAH  emission.  Besides  these  three  sources,  14  LINER  and  three 
H  ii  ULIRGs’  nuclei  show  strong  absorption  features  whose  absolute  optical  depths  suggest  an  energy  source  more 
centrally  concentrated  than  the  surrounding  dust,  such  as  a  buried  AGN.  In  total,  1 7  out  of 27  (63%)  LINER  and  3  out 
of  13  (23%)Hn  ULIRGs’  nuclei  show  some  degree  of  evidence  for  powerful  buried  AGNs,  suggesting  that  powerful 
buried  AGNs  may  be  more  common  in  LINER  ULIRGs  than  in  H  n  ULIRGs.  The  evidence  of  AGNs  is  found  in  non- 
Seyfert  ULIRGs  with  both  warm  and  cool  far-infrared  colors.  These  spectra  are  compared  with  those  of  15  ULIRGs’ 
nuclei  with  optical  Seyfert  signatures  taken  for  comparison.  The  overall  spectral  properties  suggest  that  the  total 
amount  of  dust  around  buried  AGNs  in  non-Seyfert  ULIRGs  is  systematically  larger  than  that  around  AGNs  in 
Seyfert  2  ULIRGs.  We  argue  that  the  optical  (non (detectability  of  Seyfert  signatures  in  ULIRGs  is  highly  dependent 
on  how  deeply  buried  the  AGNs  are,  and  that  it  is  essential  to  properly  evaluate  the  energetic  importance  of  buried 
AGNs  in  non-Seyfert  ULIRGs. 

Subject  headings:  galaxies:  active  —  galaxies:  ISM  —  galaxies:  nuclei  —  galaxies:  Seyfert  — 
galaxies:  starburst  —  infrared:  galaxies 


1.  INTRODUCTION 

Ultraluminous  infrared  galaxies  (  ULIRGs),  discovered  by  the 
IRAS  all  sky  survey,  radiate  quasar-like  luminosities  (>1012  L  -J) 
as  infrared  dust  emission  (Sanders  &  Mirabel  1996).  Powerful 
energy  sources,  starbursts,  and/or  active  galactic  nuclei  (AGNs) 
must  be  present  hidden  behind  dust.  This  ULIRG  population 
dominates  not  only  the  bright  end  of  the  local  luminosity  func¬ 
tion  (Soifer  et  al.  1987)  but  also  the  cosmic  submillimeter  back¬ 
ground  emission  (Blain  et  al.  1999),  so  that  their  energy  sources 
are  closely  related  to  understanding  the  connection  between  star- 
bursts  and  AGNs. 

It  is  now  generally  accepted  that  the  maj  ority  of  nearby  ULIRGs 
are  energetically  dominated  by  compact  (<500  pc),  highly  dust- 
obscured  nuclear  cores,  rather  than  extended  (>kpc),  weakly  ob¬ 
scured  starburst  activity  in  host  galaxies  (Soifer  et  al.  2000;  Fischer 
2000).  Determining  whether  the  compact  cores  are  powered  by 

1  Based  in  part  on  data  collected  at  Subaru  Telescope,  which  is  operated  by 
the  National  Astronomical  Observatory  of  Japan. 

2  Visiting  Astronomer  at  the  Infrared  Telescope  Facility,  which  is  operated  by 
the  University  of  Hawaii  under  Cooperative  Agreement  NCC  5-538  with  the 
National  Aeronautics  and  Space  Administration,  Office  of  Space  Science,  Plan¬ 
etary  Astronomy  Program. 


AGNs  and/or  by  very  compact  starbursts  is  essential,  but  it  is  not 
easy  to  distinguish  the  energy  sources  of  ULIRGs’  cores  because 
the  large  amount  of  dust  and  molecular  gas  concentrated  in  the 
nuclei  (Sanders  &  Mirabel  1996;  Spoon  et  al.  2002)  can  hide  the 
signatures  of  AGNs  and  make  their  detection  difficult.  If  AGNs 
are  obscured  by  dust  in  a  torus  geometry,  as  is  inferred  for  some 
types  of  nearby  AGN  populations  (Antonucci  1993),  clouds  along 
the  torus  axis  are  photoionized  by  the  AGN’s  hard  radiation.  At 
~  1 00  pc  distances  from  the  AGN,  these  clouds  form  the  so-called 
narrow-line  regions  (NLRs;  Robson  1996),  and  the  NLRs  pro¬ 
duce  optical  emission  lines  with  flux  ratios  that  are  distinguishable 
from  those  of  normal  starburst  galaxies  (Veilleux  &  Osterbrock 
1987).  Thus,  AGNs  obscured  by  torus-shaped  dust  distributions 
can  still  be  detectable  through  optical  spectroscopy  (classified  as 
Seyfert  2;  Veilleux  et  al.  1 999a).  Such  obscured  AGNs  with  prom¬ 
inent  NLRs  are  also  recognizable  through  strong  high-excitation 
forbidden  line  emission,  originating  in  the  NLRs,  in  high-resolution 
near-  to  mid-infrared  spectra  (Genzel  et  al.  1998;  Veilleux  et  al. 
1999b;  Murphy  et  al.  2001;  Armus  et  al.  2004).  However,  when 
virtually  all  directions  around  an  AGN  become  opaque  to  the 
bulk  of  the  AGN’s  ionizing  photons  (mostly  soft  X-rays  and  UV), 
no  significant  NLRs  develop.  In  such  buried  AGNs,  X-ray  dis¬ 
sociation  regions  (XDRs),  dominated  by  low-ionization  species 


114 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2QQ£  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2006  to  00-00-2006 

4.  TITLE  AND  SUBTITLE 

Infrared  3-4  um  Spectroscopic  Investigation  of  a  Large  Sample  of  Nearby 
Ultraluminous  Infrared  Galaxies 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory, Code  7213,4555  Overlook  Avenue, 

SW, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

unclassified  unclassified  unclassified 

24 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Z-BAND  SPECTROSCOPY  OF  ULIRGs 


115 


(Maloney  et  al.  1996),  develop,  rather  than  NLRs,  so  that  the  sig¬ 
natures  of  buried  AGNs  are  difficult  to  find  using  the  conventional 
searches  for  emission  lines  originating  in  the  NLRs. 

To  study  deeply  buried  AGNs  without  significant  NLRs,  ob¬ 
servations  at  wavelengths  of  low  dust  extinction  are  necessary. 
Infrared  Z-band  (A  =  2.8-4. 1  /rm)  spectroscopy  from  the  ground 
is  one  of  the  most  effective  ways  to  find  the  signatures  of  such 
buried  AGNs.  First,  dust  extinction  in  the  Z  band  is  substantially 
reduced  compared  to  shorter  wavelengths  (Al  ~  0.06 Ay;  Rieke 
&Lebofsky  1985;  Lutz  etal.  1996)  and  is  as  low  as  that  at  5- 13  /tm 
(Lutz  et  al.  1996).  Second,  the  contribution  from  AGN-powered 
dust  emission  to  an  observed  flux,  relative  to  stellar  emission,  in¬ 
creases  substantially  compared  with  A  <  2  /im.  Thus,  the  detec¬ 
tion  of  a  buried  AGN  becomes  more  feasible  than  at  X  <  2  /tm. 
Third,  normal  starburst  emission  and  AGN  emission  are 
clearly  distinguishable  using  Z-band  spectra  (Moorwood  1986; 
Imanishi  &  Dudley  2000).  In  a  normal  starburst  galaxy,  emission 
from  polycyclic  aromatic  hydrocarbon  (PALI)  molecules  is  ex¬ 
cited  by  far-UV  photons  in  photodissociation  regions,  the  inter¬ 
faces  between  molecular  gas  and  FI  n  regions  (Sellgren  1981; 
Normand  et  al.  1995).  Near  an  AGN,  the  PAH  molecules  can  be 
destroyed  by  strong  X-ray  radiation  from  the  AGN  (Vo it  1992; 
Siebenmorgen  et  al.  2004).  If  the  PAH  molecules  are  sufficiently 
shielded  by  obscuration  from  the  X-ray  emission  of  the  central 
AGN,  the  PAHs  can  survive.  However,  the  far-UV  emission 
from  the  AGN  will  also  be  attenuated  by  the  obscuration,  so  in  a 
pure  AGN  without  starburst  activity,  the  PAH  emission  will  be 
very  weak.  Instead,  a  featureless  continuum  from  hot  (~  1000  K), 
submicron-sized  dust  grains  (Draine  &  Lee  1984;  Mathis  & 
Whifien  1989)  heated  by  the  AGN  is  observed.  If  the  AGN  is 
deeply  buried  in  dust,  the  3.1  /tm  H20  ice  absoiption  feature 
produced  by  ice-covered  dust  and/or  the  3.4  pm  absorption 
feature  resulting  from  bare  carbonaceous  dust  should  be  ob¬ 
served  (  Imanishi  &  Dudley  2000;  Imanishi  et  al.  2001;  Imanishi 
&  Maloney  2003).  In  reality,  starbursts  and  AGNs  are  likely  to 
coexist  in  ULIRGs.  Hence,  it  is  important  to  disentangle  the  en¬ 
ergy  sources  of  such  starburst/AGN  composite  ULIRGs.  This  can 
be  done  using  observations  of  PAH  emission:  because  the  column 
densities  and  pressures  in  the  interstellar  medium  of  ULIRGs  are 
very  high,  the  starburst  regions  in  starburst/AGN  composite 
ULIRGs  should  produce  strong  PAH  emission  (Maloney  1999). 
In  a  pure  normal  starburst  galaxy,  the  equivalent  width  (=line-to- 
continuum  flux  ratio)  of  the  3.3  pm  PAH  emission  feature  should 
always  be  high,  regardless  of  dust  extinction  of  the  starburst 
(Moorwood  1986;  Imanishi  &  Dudley  2000),  because  both  the 
PAH  emission  and  continuum  emission  are  similarly  flux  atten¬ 
uated.  If  PAH  emission  is  detected  in  a  ULIRG,  but  its  equivalent 
width  is  substantially  smaller  than  in  normal  starburst  galaxies,  it 
will  strongly  suggest  that  the  PAH  equivalent  width  is  being  re¬ 
duced  by  the  contribution  to  the  observed  spectrum  of  a  PAH-free 
featureless  continuum,  produced  by  the  AGN.  Thus,  the  presence 
of  a  powerful  AGN  can  be  inferred  from  the  equivalent  width  of 
the  3.3  pm  PAH  emission  feature  and  dust  absorption  features, 
even  if  the  AGN  is  buried. 

In  addition,  using  Z-band  spectra,  we  have  another  independent 
way  to  detect  the  presence  of  powerful  buried  AGNs  in  ULIRGs: 
the  absolute  optical  depth  of  dust  absorption  features.  In  a  normal 
starburst,  the  energy  sources  (stars)  and  dust  are  spatially  well 
mixed,  while  in  a  buried  AGN,  the  energy  source  is  very  compact 
and  thus  should  be  more  centrally  concentrated  than  the  sur¬ 
rounding  dust  (Soifer  et  al.  2000;  Imanishi  &  Maloney  2003; 
Siebenmorgen  et  al.  2004).  In  the  mixed  dust/source  geometry, 
the  foreground  less  obscured  emission  (which  thus  shows  only 
weak  absorption  features)  dominates  the  observed  fluxes,  and 


therefore  there  are  upper  limits  to  the  absorption  optical  depths, 
unless  unusual  abundances  are  considered  and/or  there  is  a  very 
large  amount  of  foreground  dust  in  the  host  galaxy  (Imanishi  & 
Maloney  2003).  On  the  other  hand,  a  centrally  concentrated  en¬ 
ergy  source  can  show  arbitrarily  large  optical  depth  because  a  fore¬ 
ground  screen  dust  extinction  model  is  applicable  (  Imanishi  & 
Maloney  2003).  A  few  ULIRGs  have  been  found  to  show  sub¬ 
stantially  larger  optical  depths  than  the  upper  limits  estimated  for 
mixed  dust/source  geometry  and  thus  suggest  centrally  concen¬ 
trated  energy  sources  (e.g.,  UGC  5101  and  IRAS  08572+3915; 
Imanishi  &  Dudley  2000;  Imanishi  et  al.  200 1 ;  Imanishi  &  Maloney 
2003).  These  sources  also  display  signatures  of  buried  AGNs 
based  on  the  equivalent  width  of  the  3.3  \l m  PAH  emission 
(Imanishi  &  Dudley  2000;  Imanishi  et  al.  2001)  and  observa¬ 
tions  at  other  wavelengths  (Dudley  &  Wynn-Williams  1997; 
Soifer  et  al.  2000;  Imanishi  et  al.  2003),  demonstrating  that  this 
second  method  can  also  be  applied  to  investigate  the  energy 
sources  of  ULIRGs. 

A  number  of  papers  have  presented  Z-band  spectra  of  nearby 
ULIRGs  ( Imanishi  &  Dudley  2000;  Imanishi  et  al.  2001;  Imanishi 
&  Maloney  2003;  Risaliti  et  al.  2003,  2005),  and  their  energy 
sources  have  been  investigated.  However,  the  number  of  the  ob¬ 
served  sources  is  still  too  small  to  draw  statistically  meaningful 
conclusions.  In  this  paper  we  present  ground-based  Z-band  spec¬ 
tra  of  a  large  sample  of  nearby  ULIRGs.  Throughout  this  paper, 
H0  =  75  km  s_1  Mpc~\  =  0.3,  and  Ha  =  0.7  are  adopted. 
The  physical  scale  of  1"  is  0.34  kpc  in  the  nearest  source  at 
z  =  0.018,  0.91  kpc  at  z  =  0.05,  and  2.44  kpc  in  the  farthest 
source  atz  =  0.15. 

2.  TARGETS 

ULIRGs  in  the  IRAS  1  Jy  sample  (Kim  &  Sanders  1998),  with 
no  obvious  Seyfert  signatures  in  the  optical  (non-Seyfert  ULIRGs), 
are  our  main  targets  to  look  for  buried  AGNs.  We  restrict  the 
sample  to  those  at  z  <  0.15  so  that  the  long-wavelength  side  of 
the  redshifted  3.4  /im  dust  absorption  feature  is  at  X  <  4.1  /im; 
beyond  this  wavelength  the  steeply  increasing  Earth’s  atmo¬ 
spheric  background  signal  degrades  spectral  quality  drastically 
in  the  case  of  ground-based  observations.  A  declination  cutoff 
greater  than  —26°  was  also  applied  so  that  the  targets  are  ob¬ 
servable  with  small  air  masses  from  Mauna  Kea,  Hawaii,  our 
main  observing  site.  These  selection  criteria  result  in  24  ULIRGs 
classified  optically  as  LINERs,  1 7  ULIRGs  classified  optically 
as  H  ii  regions,  and  four  optically  unclassified  ULIRGs,  at  z  < 
0. 1 5  in  this  sample,  where  the  optical  classifications  are  based  on 
Veilleux  et  al.  (1999a;  their  Table  2). 

We  gave  highest  priority  to  LINER  ULIRGs  because  a  pure 
buried  AGN  is  expected  to  produce  XDRs  (Maloney  et  al.  1996) 
and  to  show  a  LINER-type  optical  spectrum  (Imanishi  et  al.  200 1). 
All  24  LINER  ULIRGs  were  observed,  comprising  a  complete 
sample.  For  the  H  n  and  unclassified  ULIRGs,  11  (out  of  17)  and 
two  (out  of  four)  sources  were  observed,  respectively.  One  ULIRG, 
IRAS  13443+0822,  consists  of  northeastern  and  southwestern 
nuclei,  which  are  optically  classified  as  an  H  n  region  and  a 
Seyfert  2  galaxy,  respectively  (Veilleux  et  al.  1999a).  Although 
this  ULIRG  is  classified  as  Seyfert  2  for  the  whole  system 
( Veilleux  et  al.  1 999a,  their  Table  2),  we  observed  only  the  north¬ 
eastern  H  ii  nucleus.  We  include  this  source  in  the  category  of 
H  ii  ULIRGs,  raising  the  number  of  observed  H  n  ULIRGs  to  12. 
The  observed  12  H  n  ULIRGs  also  constitute  a  complete  sample, 
in  the  right  ascension  range  10-21  hr.  The  details  of  the  observed 
non-Seyfert  ULIRGs  are  summarized  in  Table  1 . 

In  addition  to  these  non-Seyfert  ULIRGs,  several  ULIRGs 
with  obvious  Seyfert  signatures  in  the  optical  were  also  observed 


TABLE  1 

Observed  ULIRGs  and  Their  Far-Infrared  Emission  Properties 


Object 

(1) 

Redshift 

(2) 

f\2 

(Jy) 

(3) 

fl5 

(Jy) 

(4) 

,/60 

(Jy) 

(5) 

/ioo 

(Jy) 

(6) 

logiiR 

(i0) 

(7) 

fl5  If 60 
(8) 

Optical  Class 

(9) 

IRAS  00188-0856 . 

0.128 

<0.12 

0.37 

2.59 

3.40 

12.3 

0.14 

(C) 

LINER 

IRAS  03250+1606 . 

0.129 

<0.10 

<0.15 

1.38 

1.77 

12.1 

<0.11 

(C) 

LINER 

IRAS  08572+3915 . 

0.058 

0.32 

1.70 

7.43 

4.59 

12.1 

0.23 

(W) 

LINER 

IRAS  09039+0503 . 

0.125 

0.07 

0.12 

1.48 

2.06 

12.1 

0.08 

(C) 

LINER 

IRAS  09116+0334 . 

0.146 

<0.09 

<0.14 

1.09 

1.82 

12.2 

<0.13 

(C) 

LINER 

IRAS  09539+0857 . 

0.129 

<0.15 

<0.15 

1.44 

1.04 

12.0 

<0.10 

(C) 

LINER 

IRAS  10378+1108 . 

0.136 

<0.11 

0.24 

2.28 

1.82 

12.3 

0.11 

(C) 

LINER 

IRAS  10485-1447 . 

0.133 

<0.11 

0.25 

1.73 

1.66 

12.2 

0.14 

(C) 

LINER 

IRAS  10494+4424 . 

0.092 

<0.12 

0.16 

3.53 

5.41 

12.2 

0.05 

(C) 

LINER 

IRAS  11095-0238 . 

0.106 

0.06 

0.42 

3.25 

2.53 

12.2 

0.13 

(C) 

LINER 

IRAS  12112+0305 . 

0.073 

0.12 

0.51 

8.50 

9.98 

12.3 

0.06 

(C) 

LINER 

IRAS  12127-1412 . 

0.133 

<0.13 

0.24 

1.54 

1.13 

12.2 

0.16 

(C) 

LINER 

IRAS  12359-0725 . 

0.138 

0.09 

0.15 

1.33 

1.12 

12.1 

0.11 

(C) 

LINER 

IRAS  14252-1550 . 

0.149 

<0.09 

<0.23 

1.15 

1.86 

12.3 

<0.20 

(C) 

LINER 

IRAS  14348-1447 . 

0.083 

0.07 

0.49 

6.87 

7.07 

12.3 

0.07 

(C) 

LINER 

IRAS  15327+2340  (Arp  220) . 

0.018 

0.48 

7.92 

103.33 

112.40 

12.1 

0.08 

(C) 

LINER 

IRAS  16090-0139 . 

0.134 

0.09 

0.26 

3.61 

4.87 

12.5 

0.07 

(C) 

LINER 

IRAS  16468+5200 . 

0.150 

<0.06 

0.10 

1.01 

1.04 

12.1 

0.10 

(C) 

LINER 

IRAS  16487+5447 . 

0.104 

<0.07 

0.20 

2.88 

3.07 

12.2 

0.07 

(C) 

LINER 

IRAS  17028+5817 . 

0.106 

<0.06 

0.10 

2.43 

3.91 

12.1 

0.04 

(C) 

LINER 

IRAS  17044+6720 . 

0.135 

<0.07 

0.36 

1.28 

0.98 

12.1 

0.28 

(W) 

LINER 

IRAS  21329-2346 . 

0.125 

0.05 

0.12 

1.65 

2.22 

12.1 

0.07 

(C) 

LINER 

IRAS  23234+0946 . 

0.128 

<0.06 

0.08 

1.56 

2.11 

12.1 

0.05 

(C) 

LINER 

IRAS  23327+2913 . 

0.107 

<0.06 

0.22 

2.10 

2.81 

12.1 

0.10 

(C) 

LINER 

IRAS  10190+1322 . 

0.077 

<0.07 

0.38 

3.33 

5.57 

12.0 

0.11 

(C) 

H  ii 

IRAS  11387+4116 . 

0.149 

0.12 

<0.14 

1.02 

1.51 

12.2 

<0.14 

(C) 

H  ii 

IRAS  11506+1331  . 

0.127 

<0.10 

0.19 

2.58 

3.32 

12.3 

0.07 

(C) 

H  ii 

IRAS  13443+0802 . 

0.135 

<0.12 

<0.11 

1.50 

1.99 

12.2 

<0.07 

(C) 

H  na 

IRAS  13509+0442 . 

0.136 

0.10 

<0.23 

1.56 

2.53 

12.3 

<0.15 

(C) 

H  ii 

IRAS  13539+2920 . 

0.108 

<0.09 

0.12 

1.83 

2.73 

12.1 

0.07 

(C) 

H  ii 

IRAS  14060+2919 . 

0.117 

<0.10 

0.14 

1.61 

2.42 

12.1 

0.09 

(C) 

H  ii 

IRAS  15206+3342 . 

0.125 

0.08 

0.35 

1.77 

1.89 

12.2 

0.19 

(C) 

H  ii 

IRAS  15225+2350 . 

0.139 

<0.07 

0.18 

1.30 

1.48 

12.2 

0.14 

(C) 

H  ii 

IRAS  16474+3430 . 

0.111 

<0.13 

0.20 

2.27 

2.88 

12.2 

0.09 

(C) 

H  ii 

IRAS  20414-1651 . 

0.086 

<0.65 

0.35 

4.36 

5.25 

12.3 

0.08 

(C) 

H  ii 

IRAS  21208-0519 . 

0.130 

<0.09 

<0.15 

1.17 

1.66 

12.1 

<0.13 

(C) 

H  ii 

IRAS  14197+0813 . 

0.131 

<0.17 

<0.19 

1.10 

1.66 

12.2 

<0.17 

(C) 

Unclassified 

IRAS  14485-2434 . 

0.148 

<0.11 

<0.15 

1.02 

1.05 

12.2 

<0.15 

(C) 

Unclassified 

IRAS  05189-2524 . 

0.042 

0.73 

3.44 

13.67 

11.36 

12.1 

0.25 

(W) 

Sy2 

IRAS  08559+1053 . 

0.148 

<0.10 

0.19 

1.12 

1.95 

12.2 

0.17 

(C) 

Sy2 

IRAS  12072-0444 . 

0.129 

<0.12 

0.54 

2.46 

2.47 

12.4 

0.22 

(W) 

Sy2 

IRAS  13428+5608  (Mrk  273) . 

0.037 

0.24 

2.28 

21.74 

21.38 

12.1 

0.10 

(C) 

Sy2 

IRAS  13451+1232  (PKS  1345+12) . 

0.122 

0.14 

0.67 

1.92 

2.06 

12.3 

0.35 

(W) 

Sy2 

IRAS  15130-1958 . 

0.109 

<0.14 

0.39 

1.92 

2.30 

12.1 

0.20 

(W) 

Sy2 

IRAS  17179+5444 . 

0.147 

<0.08 

0.20 

1.36 

1.91 

12.3 

0.15 

(C) 

Sy2 

IRAS  01572+0009  (Mrk  1014) . 

0.163 

0.12 

0.54 

2.22 

2.16 

12.6 

0.24 

(W) 

Syl 

IRAS  07599+6508 . 

0.149 

0.26 

0.53 

1.69 

1.73 

12.5 

0.31 

(W) 

Syl 

IRAS  11598-0112 . 

0.151 

<0.14 

<0.38 

2.39 

2.63 

12.5 

<0.16 

(C) 

Syl 

IRAS  12265+0219  (3C  273) . 

0.159 

0.55 

0.90 

2.06 

2.89 

12.8 

0.44 

(W) 

Syl 

IRAS  12540+5708  (Mrk  231) . 

0.042 

1.87 

8.66 

31.99 

30.29 

12.5 

0.27 

(W) 

Syl 

IRAS  15462-0450 . 

0.100 

0.10 

0.45 

2.92 

3.00 

12.2 

0.15 

(C) 

Syl 

IRAS  21219-1757 . 

0.112 

0.21 

0.45 

1.07 

1.18 

12.1 

0.42 

(W) 

Syl 

Notes. — Col.  (1):  Object  name.  Col.  (2):  Redshift.  Cols.  (3) — (6):  f\2,f25>f60>  and/ioo  are  IRAS  fluxes  at  12, 25, 60,  and  100  fim,  respectively,  taken  from  Kim  &  Sanders 
(1998).  Col.  (7):  Decimal  logarithm  of  infrared  (8-1000  /im)  luminosity  in  units  of  solar  luminosity  (Z,Q),  calculated  with  Ljr  =  (2.1  x  1039)Z)(Mpc)2(13.48/i2  + 
5. 16/25  +  2.58^60  +  /ioo)  ergs  s-1  (Sanders  &  Mirabel  1996).  Since  the  calculation  is  based  on  our  adopted  cosmology,  the  infrared  luminosities  slightly  (<10%)  differ 
from  the  values  shown  in  Kim  &  Sanders  (1998,  their  Table  1 ,  col.  [  1 5]).  For  sources  that  have  upper  limits  in  some  IRAS  bands,  we  assume  that  the  actual  flux  is  the  upper 
limit.  Even  though  zero  values  are  adopted,  the  difference  of  the  infrared  luminosities  is  very  small,  less  than  0.2  dex,  which  will  not  affect  our  main  conclusions  in  this 
paper.  Col.  (8):  ULIRGs  with/^/^o  <  0-2  and  >0.2  are  classified  as  cool  and  warm  (denoted  as  “C”  and  “ W”),  respectively  (Sanders  et  al.  1 988).  Col.  (9):  Optical  spectral 
classification  by  Veilleux  et  al.  (1999a). 

a  The  northeastern  nucleus  (Kim  et  al.  2002)  classified  optically  as  an  H  n  region  (Veilleux  et  al.  1999a)  is  observed. 


Z-BAND  SPECTROSCOPY  OF  ULIRGs 


117 


as  a  control  sample.  There  are  12  Seyfert  2  ULIRGs  atz  <  0.15 
and  declination  greater  than  —26°  in  the  IRAS  1  Jy  sample, 
excluding  the  above-mentioned  IRAS  13443+0822.  In  total,  7 
Seyfert  2  ULIRGs  (Table  1)  were  observed.  The  main  aim  of 
observing  these  Seyfert  2  ULIRGs  is  twofold:  (1)  to  check  if  our 
Z-band  spectroscopic  method  can  successfully  find  signatures  of 
obscured  AGNs  in  these  ULIRGs  that  we  know  possess  obscured 
AGNs,  and  (2)  to  provide  a  control  sample  with  which  the  prop¬ 
erties  of  buried  AGNs  in  non-Seyfert  ULIRGs  can  be  compared. 

Z-band  spectra  of  several  Seyfert  1  ULIRGs  were  also  ob¬ 
tained.  For  these  unobscured  AGNs,  since  no  significant  3.4  pm 
dust  absorption  feature  is  expected,  we  extend  the  redshift  restric¬ 
tion  to  z  <  0. 1 7  where  the  longer  wavelength  side  of  the  3.3  pm 
PALI  emission  feature  is  covered  in  the  ground-based  Z-band 
spectrum.  There  are  7  Seyfert  1  ULIRGs  atz  <  0.17  with  decli¬ 
nation  greater  than  —26°  in  the  IRAS  1  Jy  sample,  and  all  of  these 
sources  were  observed  (Table  1). 

3.  OBSERVATIONS  AND  DATA  ANALYSIS 

Observations  were  made  primarily  using  the  1RCS  near-infrared 
spectrograph  (Kobayashi  et  al.  2000)  on  the  Subaru  8.2  m  tele¬ 
scope  (lye  et  al.  2004).  Several  bright  ULIRGs  were  observed  with 
CGS4  (Mountain  et  al.  1990)  on  the  UKIRT  3.6  m  telescope  and 
with  SpeX  (Rayner  et  al.  2003)  and  NSFCAM  (Shure  et  al.  1994) 
on  the  Infrared  Telescope  Facility  (IRTF)  3  m  telescope.  Table  2 
provides  a  detailed  observing  log. 

For  the  Subaru  1RCS  observing  runs,  the  sky  was  clear  during 
the  observations  of  all  but  one  UL1RG  (IRAS  16090—0139). 
The  seeing  at  K,  measured  in  images  taken  before  Z-band 
spectroscopy,  was  0"4-0"9  FWHM.  A  Of 9  wide  slit  and  the  Z 
grism  were  used  with  a  58  mas  pixel  scale.  The  achievable  spec¬ 
tral  resolution  is  R  ~  140  at  1  ~  3.5  pm.  A  standard  telescope 
nodding  technique  (ABBA  pattern),  with  a  throw  of  7"- 10" 
along  the  slit,  was  employed  to  subtract  background  emission. 
The  optical  guider  of  the  Subaru  telescope  was  used  to  monitor 
the  telescope  tracking.  Exposure  time  was  1.5-2. 5  s,  and  15-40 
co-adds  were  made  at  each  nod  position. 

The  spectra  of  two  bright  ULIRGs,  IRAS  08572+3915  and 
Mrk  1014,  were  taken  using  IRTF  SpeX  (Rayner  et  al.  2003). 
The  1. 9-4.2  pm  cross-dispersed  mode  with  a  1"6  wide  slit  was 
employed.  This  mode  enables  Z-  (2.8-4. 1  pm)  and  A+band  (2— 
2.5  pm)  spectra  to  be  obtained  simultaneously,  with  a  spectral 
resolution  of  R  ~  500.  The  sky  conditions  were  photometric 
throughout  the  observations,  and  the  seeing  at  K  was  measured  to 
be  in  the  range  0"5-0"8  FWHM.  A  standard  telescope  nodding 
technique  (ABBA  pattern)  with  a  throw  of  7"5  was  employed 
along  the  slit.  The  telescope  tracking  was  monitored  with  the 
infrared  slit  viewer  of  SpeX.  Each  exposure  was  15  s,  and  two 
co-adds  were  made  at  each  position. 

For  the  IRTF  NSFCAM  observations  oflRAS  07599+6508, 3C 
273,  and  IRAS  15462-0450,  we  used  the  NSFCAM  grism  mode 
(Shure  et  al.  1994).  Sky  conditions  were  photometric  throughout 
the  observations,  and  the  seeing  FWHM  was  measured  to  be  0"5- 
0"9  in  the  K  band.  The  HKL  grism  and  Z  blocker  were  used  with 
the  4  pixel  slit  (=1"2).  The  resulting  spectral  resolution  was  R  ~ 
150  at  X  ~  3.5  pm.  A  standard  10"  telescope  nodding  technique 
(ABBA  pattern)  was  employed.  Each  exposure  was  1  s,  and  30  co¬ 
adds  were  made  at  each  nod  position.  An  optical  guide  star  was 
used  for  the  telescope  tracking,  and  the  objects  were  observed  in  the 
eastern  sky,  where  the  tracking  performance  was  well  confirmed. 

Spectra  of  four  ULIRGs  (Arp  220,  IRAS  05189-2524,  Mrk 
273,  Mrk  231)  were  taken  with  UKIRT  CGS4  (Mountain  et  al. 
1990),  using  a  1"2  wide  slit,  and  were  published  in  Imanishi  & 
Dudley  (2000).  The  observing  details  are  found  in  that  paper. 


For  all  the  observing  runs,  A-,  F-,  and  G-type  main-sequence 
stars  (Table  2)  were  observed  as  standard  stars,  with  mean  air- 
mass  difference  of  <0.1  to  the  individual  ULIRGs’  nuclei,  to  cor¬ 
rect  for  the  transmission  of  the  Earth’s  atmosphere  and  to  provide 
flux  calibration.  The  Z-band  magnitudes  of  the  standard  stars  were 
estimated  from  their  Uband  (A  =  0.6  pm)  magnitudes,  adopting 
the  V  —  L  colors  appropriate  to  the  stellar  types  of  individual 
standard  stars  (Tokunaga  2000). 

Standard  data  analysis  procedures  were  employed,  using 
IRAF.3  Initially,  frames  taken  with  an  A  (or  B)  beam  were  sub¬ 
tracted  from  frames  subsequently  taken  with  a  B  (or  A)  beam, 
and  the  resulting  subtracted  frames  were  added  and  divided  by  a 
spectroscopic  flat  image.  Then,  bad  pixels  and  pixels  hit  by  cos¬ 
mic  rays  were  replaced  with  the  interpolated  values  of  the  sur¬ 
rounding  pixels.  Finally,  the  spectra  of  ULIRGs’  nuclei  and 
standard  stars  were  extracted,  by  integrating  signals  over  0"9- 
2"0,  depending  on  actual  signal  profiles.  Wavelength  calibration 
was  performed  taking  into  account  the  wavelength-dependent 
transmission  of  the  Earth’s  atmosphere.  The  spectra  of  ULIRGs’ 
nuclei  were  divided  by  the  observed  spectra  of  standard  stars, 
multiplied  by  the  spectra  of  blackbodies  with  temperatures  ap¬ 
propriate  to  individual  standard  stars  (Table  2). 

Flux  calibration  for  all  ULIRGs  but  one  (IRAS  16090—0139) 
was  done  based  on  signals  of  ULIRGs  and  standard  stars  de¬ 
tected  inside  our  slit  spectra.  For  the  Subaru  1RCS  spectrum  of 
IRAS  16090—0139,  since  its  spectrum  was  taken  under  non¬ 
photometric  conditions,  we  calibrated  its  flux  level  based  on  our 
own  Z-band  photometry  (Z  =  13.3),  made  on  a  different  night 
using  IRTF  NSFCAM.  Seeing  sizes  at  K  (and  Z)  were  always 
smaller  than  the  employed  slit  widths,  and  good  telescope  track¬ 
ing  performances  of  Subaru,  IRTF,  and  UKIRT  were  established. 
We  thus  expect  possible  slit  loss  to  be  minimal.  To  estimate  this 
ambiguity,  we  divided  the  whole  data  set  into  subgroups  and 
compared  their  flux  levels.  With  the  exception  of  the  highly 
time  variable  source  3C  273,  the  flux  levels  agree  within  10% 
for  most  objects  and  30%  for  the  worst  case,  suggesting  that  the 
flux  uncertainty  is  at  a  similar  level,  which  will  not  affect  our 
main  conclusions  significantly.  For  3  C  273,  time  variation  was 
clearly  recognizable  even  during  our  IRTF  NSFCAM  observing 
run.  For  this  source,  the  resulting  spectrum  is  a  time-averaged 
one. 

To  obtain  an  adequate  signal-to-noise  ratio  in  each  element, 
appropriate  binning  of  spectral  elements  was  performed,  partic¬ 
ularly  at  A0bs  <3.3  pm  in  the  observed  frame,  where  the  scatter 
of  data  points  is  larger,  due  to  poorer  Earth  atmospheric  trans¬ 
mission,  than  at  A0bS  >3.4  pm.  The  resulting  spectral  resolution 
at  A0bs  <3.3  pm  is  R Si  100  in  most  cases.  The  atmospheric 
transmission  curve  at  A0bS  =  2. 8-3. 3  pm  is  fairly  smooth  at  this 
spectral  resolution.  Thus,  even  if  the  net  positions  of  the  target 
object  and  standard  star  on  the  slit  differ  slightly  (on  the  subpixel 
scale)  along  the  wavelength  direction,  the  standard  data  analy¬ 
sis  described  above  is  expected  to  produce  no  significant  spuri¬ 
ous  features  in  spectra  with  R  Si  100  (see  Imanishi  &  Maloney 
2003). 

4.  RESULTS 
4.1.  L-Band  Spectra 

Figures  1,  2,  3,  4,  and  5  present  flux-calibrated  Z-band 
spectra  of  the  nuclei  of  ULIRGs  optically  classified  as  LINERs, 

3  IRAF  is  distributed  by  the  National  Optical  Astronomy  Observatory,  which  is 
operated  by  the  Association  of  Universities  for  Research  in  Astronomy  (AURA), 
Inc.,  under  cooperative  agreement  with  the  National  Science  Foundation. 


TABLE  2 
Observing  Log 


Standard  Stars 


Object 

a) 

Date 

(UT) 

(2) 

Telescope  Instrument 

(3) 

Integration 

(minutes) 

(4) 

P.A.a 

(deg) 

(5) 

Name 

(6) 

L  Magnitude 

(7) 

Type 

(8) 

Teff 

(K) 

(9) 

IRAS 

00188-0856 . 

2002  Aug  19 

Subaru  IRCS 

25 

90 

HR  72 

5.0 

GO  V 

5930 

IRAS 

03250+1606 . 

2002  Aug  19 

Subaru  IRCS 

23 

0 

HR  763 

4.3 

F7  V 

6240 

IRAS 

08572+39 15NWb . 

2003  Mar  18 

IRTF  SpeX 

30 

0 

HR  3451 

5.0 

F7  V 

6240 

IRAS 

09039+0503 . 

2005  Feb  21 

Subaru  IRCS 

48 

0 

HR  3492 

4.4 

A0  V 

9480 

IRAS 

091 16+03  34W . 

2005  Feb  20 

Subaru  IRCS 

40 

0 

HR  3651 

6.1 

A0  V 

9480 

IRAS 

09539+0857 . 

2005  Feb  22 

Subaru  IRCS 

64 

0 

HR  3651 

6.1 

A0  V 

9480 

IRAS 

10378+1108 . 

2005  Feb  21 

Subaru  IRCS 

32 

0 

HR  4079 

5.3 

F6  V 

6400 

IRAS 

10485-1447 . 

2005  Feb  21 

Subaru  IRCS 

32 

0 

HR  4158 

4.3 

F7  V 

6240 

IRAS 

10494+4424 . 

2005  Feb  20 

Subaru  IRCS 

32 

0 

HR  4098 

5.0 

F9  V 

6000 

IRAS 

11095-0238 . 

2005  Feb  21 

Subaru  IRCS 

32 

40c 

HR  4356 

5.4 

A0  V 

9480 

IRAS 

12112+03  05NE/SW . 

2004  Feb  8 

Subaru  IRCS 

32 

38d 

HR  4533 

4.8 

F7  V 

6240 

IRAS 

12127— 1412NE . 

2005  Feb  21 

Subaru  IRCS 

16 

0e 

HR  4529 

4.9 

F7  V 

6240 

IRAS 

12359— 0725N . 

2005  Feb  23 

Subaru  IRCS 

40 

0f 

HR  4533 

4.8 

F7  V 

6240 

IRAS 

14252- 1550E/W . 

2005  Feb  22 

Subaru  IRCS 

80 

125g 

HR  5355 

5.9 

A0 

9480 

IRAS 

14348- 144 7NE/SW . 

2005  May  28 

Subaru  IRCS 

32 

30h 

HR  5652 

4.5 

A0  V 

9480 

IRAS 

15327+2340  (Arp  220)' . 

2000  Feb  20 

UKIRT  CGS4 

64 

90 

HR  5634 

3.8 

F5  V 

6500 

IRAS 

16090-0139 . 

2003  May  18 

Subaru  IRCS 

14 

0 

HR  5859 

5.6 

A0  V 

9480 

IRAS 

16468+5200E  . 

2005  Feb  23 

Subaru  IRCS 

72 

85^ 

HR  5949 

6.3 

A0  V 

9480 

IRAS 

16487+5447W . 

2002  Mar  28 

Subaru  IRCS 

36 

70k 

HR  5949 

6.3 

A0  V 

9480 

IRAS 

17028+5817W . 

2004  Jun  12 

Subaru  IRCS 

20 

175* 

HR  5949 

6.3 

A0  V 

9480 

IRAS 

17044+6720 . 

2005  Jun  13 

Subaru  IRCS 

26 

90 

HR  6025 

5.4 

AO  V 

9480 

IRAS 

21329-2346 . 

2002  Oct  24 

Subaru  IRCS 

32 

90 

HR  7898 

4.5 

G8  V 

5400 

IRAS 

23234+0946W . 

2002  Oct  24 

Subaru  IRCS 

32 

120m 

HR  8653 

4.6 

G8  IV 

5400 

IRAS 

23327+2913S . 

2002  Oct  24 

Subaru  IRCS 

24 

90" 

HR  8955 

5.1 

F6  V 

6400 

IRAS 

10190+1322E/W . 

2005  Feb  22 

Subaru  IRCS 

32 

65° 

HR  3998 

5.1 

F7  V 

6240 

IRAS 

11387+4116 . 

2005  Feb  23 

Subaru  IRCS 

32 

0 

HR  4345 

5.0 

GO  V 

5930 

IRAS 

11506+1331  . 

2005  Feb  22 

Subaru  IRCS 

32 

0 

HR  4437 

4.7 

GO  V 

5930 

IRAS 

13443+0802NE . 

2005  May  28 

Subaru  IRCS 

24 

170p 

HR  5243 

4.9 

F6  V 

6400 

IRAS 

13509+0442 . 

2005  Feb  22 

Subaru  IRCS 

32 

0 

HR  5011 

3.8 

GO  V 

5930 

IRAS 

13539+2920NW . 

2004  Jun  12 

Subaru  IRCS 

20 

0q 

HR  5346 

4.8 

F8  V 

6000 

IRAS 

14060+2919 . 

2005  Jun  12 

Subaru  IRCS 

18 

0 

HR  5346 

4.8 

F8  V 

6000 

IRAS 

15206+3342 . 

2005  Jun  13 

Subaru  IRCS 

18 

0 

HR  5728 

4.5 

G3  V 

5800 

IRAS 

15225+2350 . 

2005  Feb  20 

Subaru  IRCS 

32 

0 

HR  5630 

5.0 

F8  V 

6000 

IRAS 

16474+3430S . 

2005  Feb  21 

Subaru  IRCS 

32 

1 6  lr 

HR  6064 

5.2 

G1  V 

5900 

IRAS 

20414-1651 . 

2005  May  28 

Subaru  IRCS 

24 

0 

HR  7855 

4.9 

F6  V 

6400 

IRAS 

21208— 0519N . 

2005  May  28 

Subaru  IRCS 

20 

17s 

HR  8041 

4.7 

G1  V 

5900 

IRAS 

14197+0813 . 

2005  May  28 

Subaru  IRCS 

32 

0 

HR  5307 

5.1 

F7  V 

6240 

IRAS 

14485-2434 . 

2005  May  28 

Subaru  IRCS 

24 

0 

HR  5504 

5.0 

F7  V 

6240 

IRAS 

05189-2524' . 

1999  Sep  9 

UKIRT  CGS4 

53 

19 

HR  1762 

4.7 

AO  V 

9480 

IRAS 

08559+1053 . 

2005  Feb  23 

Subaru  IRCS 

24 

0 

HR  3510 

4.9 

G1  V 

5900 

IRAS 

12072— 0444N/S . 

2005  May  29 

Subaru  IRCS 

16 

O' 

HR  4533 

4.8 

F7  V 

6240 

IRAS 

13428+5608  (Mrk  273)' . 

2000  Feb  20 

UKIRT  CGS4 

43 

35 

HR  4761 

4.8 

F6-8  V 

6200 

IRAS 

13451+1232W  (PKS  1345+12) . 

2005  Feb  23 

Subaru  IRCS 

32 

105“ 

HR  5011 

3.8 

GO  V 

5930 

IRAS 

15130-1958 . 

2005  May  28 

Subaru  IRCS 

12 

0 

HR  5504 

5.0 

F7  V 

6240 

IRAS 

17179+5444 . 

2005  May  28 

Subaru  IRCS 

16 

0 

HR  6360 

4.5 

G5  V 

5700 

IRAS 

01572+0009  (Mrk  1014) . 

2003  Sep  8 

IRTF  SpeX 

160 

0 

HR  650 

4.1 

F8  V 

6000 

IRAS 

07599+6508 . 

2004  Apr  4 

IRTF  NSFCAM 

40 

0 

HR  3028 

4.8 

F6  V 

6400 
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TABLE  2 — Continued 


Standard  Stars 


Date  Integration  P.A.a 

Object  (UT)  Telescope  Instrument  (minutes)  (deg)  Name  Z  Magnitude  Type  (K) 

(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9) 


IRAS  11598-0112 .  2005  May  29  Subaru  IRC S  8  0  HR  4533  4.8  F7  V  6240 

IRAS  12265+0219  (3C  273) .  2004  Apr  5  IRTF  NSFCAM  28  0  HR  4708  5.0  F8  V  6000 

IRAS  12540+5708  (Mrk  231  )v .  2000  Feb  20  UKIRT  CGS4  21  0  HR  4761  4.8  F6-8  V  6200 

IRAS  15462-0450 .  2004  Apr  5  IRTF  NSFCAM  120  0  HR  5779  5.2  F7  V  6240 

IRAS  21219-1757 .  2005  May  29  Subaru  IRC S  8  0  HR  7994  4.9  G1  V  5900 


Notes. — Col.  (1):  Object  name.  Col.  (2):  Observing  date  in  UT.  Col.  (3):  Telescope  and  instrument.  Col.  (4):  Net  on-source  integration  time  in  minutes.  Col.  (5): 
Position  angle  of  the  slit.  Col.  (6):  Standard  star  name.  Col.  (7):  Adopted  Z-band  magnitude.  Col.  (8):  Stellar  spectral  type.  Col.  (9):  Effective  temperature. 

a  0°  corresponds  to  the  north-south  direction.  Position  angle  increases  counterclockwise  on  the  sky  plane. 

b  A  spectrum  at  X  =  3.2— 3.8  fim  was  shown  by  Imanishi  &  Dudley  (2000).  We  present  a  new  spectrum  with  a  wider  wavelength  coverage  of  2  =  2.8— 4.1  fim,  to 
investigate  a  broad  3.1  /im  H2O  ice  absorption  feature.  We  observed  the  northwestern  nucleus,  which  is  much  brighter  in  the  infrared  Z  band  than  the  southeastern 
nucleus  (Zhou  et  al.  1993). 

c  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~0"5  (Bushouse  et  al.  2002).  In  our  spectrum,  a  double  peak  along  the  slit  direction  is 
recognizable,  but  signals  from  the  double  nuclei  are  not  clearly  resolvable.  A  spectrum  combining  both  nuclei  is  presented  in  this  paper. 

d  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~3"  (Scoville  et  al.  2000;  Kim  et  al.  2002). 

e  Only  the  brighter  northeastern  nucleus  (Kim  et  al.  2002)  was  observed. 

f  Only  the  northern  primary  nucleus  (Dinh-V-Trung  et  al.  2001)  was  observed. 

8  The  southwestern  nucleus  (Kim  et  al.  2002)  was  resolved  into  double  nuclei  (eastern  and  western)  with  a  separation  of  ~l''l  in  our  Z^-band  image.  Both  the 
eastern  and  western  nuclei  were  observed. 

h  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~3"  (Scoville  et  al.  2000;  Kim  et  al.  2002;  Charmandaris  et  al.  2002). 

1  The  same  spectrum  as  previously  shown  by  Imanishi  &  Dudley  (2000)  is  presented  again.  The  wavelength  coverage  is  narrower  than  the  other  ULIRGs. 

J  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~3"5  (Kim  et  al.  2002).  Only  the  eastern  nucleus  was  bright  enough  to  extract  a  meaningful 
Z-band  spectrum. 

k  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~3"  (Murphy  et  al.  1996;  Kim  et  al.  2002).  Only  the  western  primary  nucleus  (Murphy  et  al. 
2001)  was  bright  enough  to  extract  a  meaningful  Z-band  spectrum. 

1  Only  the  western  primary  nucleus  (Dinh-V-Trung  et  al.  2001)  was  observed. 

m  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~4"  in  the  ZT-band  image  (Kim  et  al.  2002).  Only  the  western  nucleus  was  bright  enough  to 
extract  a  meaningful  Z-band  spectrum. 

n  Only  the  southern  primary  nucleus  (Dinh-V-Trung  et  al.  2001)  was  observed. 

0  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~4"  (Kim  et  al.  2002). 

p  To  simultaneously  observe  northeastern  and  eastern  nuclei  with  a  separation  of  ~5"  (Kim  et  al.  2002).  Both  nuclei  were  detected,  but  only  the  northeastern 
nucleus  was  bright  enough  to  extract  a  meaningful  Z-band  spectrum. 

q  Only  the  northwestern  primary  nucleus  (Kim  et  al.  2002)  was  observed. 

r  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~3?5  (Kim  et  al.  2002).  Only  the  southern  nucleus  was  bright  enough  to  extract  a  meaningful 
Z-band  spectrum. 

s  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~7"  (Kim  et  al.  2002).  Only  the  northern  nucleus  was  bright  enough  to  extract  a  meaningful 
Z-band  spectrum. 

1  The  nucleus  (Kim  et  al.  2002)  was  resolved  into  double  nuclei  (northern  and  southern)  with  a  separation  of  ~0"9  in  our  Z>band  image.  Both  nuclei  were 
observed. 

u  To  simultaneously  observe  double  nuclei  with  a  separation  of  ~2"  (Scoville  et  al.  2000).  Only  the  western  AGN-harboring  main  nucleus  (Imanishi  & 
Terashima  2004)  was  bright  enough  to  extract  a  meaningful  Z-band  spectrum. 

v  Z-band  spectra  have  previously  been  presented  by  Imanishi  et  al.  (1998)  and  Imanishi  &  Dudley  (2000).  The  spectrum  by  Imanishi  &  Dudley  (2000)  is  shown 
again.  The  wavelength  coverage  is  narrower  than  the  other  ULIRGs. 


H 11  regions,  unclassified,  Seyfert  2,  and  Seyfert  1,  respectively. 
For  ULIRGs  that  have  more  than  one  nuclei  brighter  than  Z  ~ 
14.5  mag,  spectra  of  individual  nuclei  are  shown  separately.  The 
total  number  of  Z-band  spectra  is  27  for  LINER,  1 3  for  H  n,  two 
for  unclassified,  eight  for  Seyfert  2,  and  seven  for  Seyfert  1 
ULIRGs’  nuclei.  Z-band  spectra  of  several  ULIRGs’  nuclei,  pre¬ 
viously  presented  by  Imanishi  &  Dudley  (2000)  and  Imanishi  & 
Maloney  (2003),  are  shown  again  here. 

There  are  a  variety  of  spectral  shapes,  but  the  majority  of 
the  observed  ULIRGs  show  emission  features  at  20bS  =  (1  +  z)  x 
(3.29  /mi),  the  wavelength  where  the  3.3  fim  PAH  emission 
feature  has  a  peak.  We  thus  identify  these  features  as  the  3.3  fim 
PAH  emission.  To  estimate  the  strength  of  the  3.3  fim  PAH 
emission  feature  in  the  spectra  where  broad  absorption  features 
coexist  and  overlap  with  the  PAH  emission  feature  in  wave¬ 
length,  we  make  the  reasonable  assumption  that  the  profiles  of 
the  3.3  /imPAH  emission  in  these  ULIRGs  are  similar  to  those  of 
Galactic  star-forming  regions  and  nearby  starburst  galaxies;  the 
main  emission  profile  extends  between  Arest  =  3.24  and  3.35  fan 


in  the  rest  frame  (Tokunaga  et  al.  1991;  Imanishi  &  Dudley 
2000).  Data  points  at  slightly  shorter  than  2rest  =  3.24  fim  and 
slightly  longer  than  /lrest  =  3.35  fim,  unaffected  by  obvious  ab¬ 
sorption  features,  are  adopted  as  the  continuum  levels  to  estimate 
the  3.3  fim  PAH  emission  strength.  We  adopt  the  spectral  profile 
of  type  1  sources  (Tokunaga  et  al.  1991)  as  a  template  for  the 
3.3  fim  PAH  emission.  The  adopted  template  reproduces  the  ob¬ 
served  3.3  fim  PAH  emission  features  of  the  ULIRGs  reasonably 
well,  with  our  spectral  resolution  and  signal-to-noise  ratios.  Table  3 
summarizes  the  fluxes,  luminosities,  and  rest-frame  equivalent 
widths  (EW3  3pah)  of  the  3.3  fim  PAH  emission  feature.  The  uncer¬ 
tainties  of  the  3 .3  fim  PAH  emission  fluxes,  estimated  from  the  fit¬ 
tings,  are  <20%  in  most  cases  and  ~30%  even  for  ULIRGs  with 
large  scatters  at  the  wavelengths  around  the  3.3  fim  PAH  emis¬ 
sion.  We  also  tried  several  continuum  levels  in  a  reasonable  range, 
to  estimate  possible  systematic  uncertainties  coming  from  contin¬ 
uum  determination  ambiguities.  The  estimated  3.3  fim  PAH  emis¬ 
sion  fluxes  usually  agree  within  20%  for  ULIRGs  with  smooth 
continuum  emission.  Even  in  ULIRGs  with  winding  spectra,  the 
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Fig.  1. — Infrared  L-band  (lobs  —  2.8— 4.1  fim)  spectra  of  ULIRGs’  nuclei  classified  optically  as  LINERs.  The  abscissa  and  ordinate  are  the  observed  wavelength  in 
fim  and  Fx  in  10~16  W  m-2  /tm_1,  respectively.  The  lower  arrows  with  “3.3  PAH”  indicate  the  expected  wavelength  of  the  3.3  fim  PAH  emission  (2rest  =  3.29  fin i). 
The  solid  lines  are  the  adopted  continuum  to  estimate  the  3.3  fim  PAH  emission  fluxes.  The  arrows  with  “3.05  /im”  and  “3.4  fin i”  indicate  the  expected  absorption  peak 
at  Irest  —  3.05  and  3.4  fim,  respectively.  For  ULIRGs  with  obvious  signatures  of  these  absorption  features,  dashed  straight  lines  are  plotted  to  show  the  adopted 
continuum  levels  used  to  estimate  the  absorption  optical  depths.  The  dotted  lines  indicate  the  wavelength  range  where  the  effects  of  the  broad  3.05  fim  ice  absorption 
feature  can  be  significant  (2rest  =  2.75—3.55  fim  adopted  from  the  spectrum  of  Elias  16;  Smith  et  al.  1989).  The  ice  absorption  feature  is  usually  very  strong  at 
lrest  =  2.9— 3.2  fim,  but  the  profiles  and  wavelength  range  of  weaker  absorption  wings  are  found  to  vary  among  different  Galactic  objects  (Smith  et  al.  1989). 


estimated  fluxes  do  not  differ  by  >30%,  as  long  as  reasonable  con¬ 
tinuum  levels  are  adopted.  Thus,  the  total  PAH  flux  uncertainties 
are  unlikely  to  exceed  ^40%  even  in  the  worst  case. 

In  addition  to  the  3 .3  fim  PAH  emission  feature,  clear  absorption 
features  at  Arest  =  3.4  fim  by  bare  carbonaceous  dust  (Pendleton 
et  al.  1994)  were  found  in  three  LINER  ULIRGs,  IRAS  08572+ 
3 9 1 5NW,  IRAS  1 2 1 27- 1 4 1 2NE,  and  IRAS  1 7044+6720  ( Fig.  1 ), 
and  one  Seyfert  2  UL1RG,  IRAS  12072-0444N  (Fig.  4).  The 
unobscured  Z-band  spectral  shapes  of  AGNs  are  dominated  by  hot 
dust  emission  and  almost  linear  (Imanishi  2002, 2003;  Imanishi  & 


Wada  2004).  Starburst  galaxies  also  show  Z-band  spectra  that  are 
approximately  linear,  aside  from  the  3.3  fim  PAH  emission  feature 
(Imanishi  &  Dudley  2000).  To  estimate  the  optical  depths  of  the 
3.4  fim  carbonaceous  dust  absorption  features,  we  therefore  adopt 
a  linear  continuum,  which  is  shown  by  a  dashed  line  in  the  spectra 
of  these  four  ULIRGs  in  Figures  1  and  4.  The  estimated  optical 
depths  of  the  3.4  fim  absorption  feature  (7-3.4)  for  these  ULIRGs 
are  summarized  in  Table  4. 

In  the  observed  spectra  of  many  ULIRGs,  other  important 
features  are  present:  ( 1 )  There  is  a  spectral  “gap,”  in  that  the  flux 
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Fig.  1. — Continued 


densities  at  the  shorter  wavelength  side  of  the  3.3  pm  PAH  emis¬ 
sion  feature  are  fainter  than  the  extrapolation  from  the  data  points 
at  the  longer  wavelength  side.  (2)  The  continuum  is  concave.  At 
the  shorter  wavelength  side  of  the  3.3  /mi  PAH  emission,  the 
continuum  flux  level  initially  decreases  with  decreasing  wave¬ 
length  but  then  begins  to  increase  again  at  lrest  ~  3.05  pm,  and 
then  the  spectra  become  flat  at  the  shortest  wavelengths  (Arest  < 
2.7  pm).  This  behavior  is  naturally  explained  by  the  broad  H20 
ice  absorption  feature  caused  by  ice-covered  dust  grains  (Spoon 
et  al.  2000;  Imanishi  &  Maloney  2003).  The  signatures  of  this 
feature  are  found  in  none  of  the  Seyfert  1  ULIRGs,  as  expected 
from  the  unobscured  view  of  the  AGN.  The  emission  from  stars 
is  reduced  in  the  Z  band  compared  to  that  in  the  K  band,  whereas 
hot  (^1000  K)  dust  heated  by  an  AGN  has  strong  emission  at 
both  K  and  Z.  In  a  starburst/AGN  composite  ULIRG,  if  the 
Z-band  emission  were  dominated  by  the  AGN-powered  hot  dust 


emission,  while  the  longer  wavelength  tail  of  the  stellar  emission 
extended  to  the  shorter  part  of  the  Z-band  spectrum,  the  Z-band 
spectrum  may  appear  to  be  concave.  However,  this  scenario  can¬ 
not  produce  the  observed  spectral  gaps  (point  1  above),  so  the 
detected  features  are  ascribed  to  H20  ice  absorption. 

The  shortest  wavelength  of  the  H20  ice  absorption  feature  is 
4rest  ~  2.7  pm  ( Whittet  et  al.  1988;  Smith  et  al.  1993).  Due  to 
redshifting,  for  the  majority  of  the  observed  ULIRGs  at  z  >  0.07, 
the  continuum  on  the  short-wavelength  side  of  the  absorption 
feature  is  included  in  our  Z-band  spectra  at  /lobs  ~  2.8-4. 1  pm. 
Thus,  continuum  levels  at  both  shorter  and  longer  wavelength 
sides  of  important  features  (broad  3.1  pm  ice  absorption,  3.3  pm 
PAH  emission,  and  3.4  pm  absorption)  are  covered  in  our  spec¬ 
tra,  making  continuum  determination  ambiguities  small.  To 
estimate  the  optical  depth  of  this  H20  absorption  feature,  we 
adopt  a  linear  continuum,  as  we  did  for  the  3.4  pm  absorption 
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Fig.  1. — Continued 


feature.  Although  Imanishi  &  Maloney  (2003)  adopted  a  con¬ 
cave  quadratic  continuum  to  derive  the  lower  limits,  linear 
continua  are  arguably  more  plausible,  given  that  both  the  AGN 
and  starburst  components  (outside  the  3.3  pm  PAH  emission 
feature)  show  linear  continuum  emission  (Imanishi  &  Dudley 
2000;  Imanishi  2002,  2003;  Imanishi  &  Wada  2004).  This  H20 
ice  absorption  feature  is  spectrally  very  broad,  and  the  con¬ 
tinuum  level  can  vary  slightly,  depending  on  the  adopted  data 
points  used  for  the  continuum  determination.  We  tried  several 
plausible  linear  continuum  levels  and  finally  adopt  the  lowest 
plausible  level  to  obtain  a  conservative  lower  limit  for  the  op¬ 
tical  depth  of  the  3.1  //m  H20  ice  absorption  feature  (T3 1).  The 
derived  T3  1  values  for  ULIRGs  showing  this  feature  clearly  are 
summarized  in  Table  4. 


4.2.  K-Band  Spectra 

Two  ULIRGs,  IRAS  08572+3915  and  Mrk  1014,  were  ob¬ 
served  with  1RTF  SpeX,  and  so  A'- band  (A  =  2. 0-2. 5  //m)  spec¬ 
tra  were  simultaneously  taken.  Figure  6  shows  the  X-band  spectra 
of  these  two  sources. 

The  A'-band  spectrum  oflRAS  08572+3915  (Fig.  6,  left  panel) 
is  similar  to  that  presented  by  Goldader  et  al.  (1995).  The  spectrum 
is  very  red  and  displays  no  detectable  CO  absorption  features  due 
to  stars  at  Arest  >  2.29  pm  or  20bS  >  2.43  pm,  which  suggests  that 
featureless  continuum  emission  from  dust  heated  by  an  AGN  dom¬ 
inates  the  A'-band  spectrum. 

The  X-band  spectrum  of  Mrk  1014  (Fig.  6,  right  panel)  shows 
a  strong  Paa  emission  line.  The  emission  line  displays  both 


No.  1,  2006 


I-BAND  SPECTROSCOPY  OF  ULIRGs 


123 


observed  wavelength  (/im) 


Fig.  1. — Continued 


broad  and  narrow  components,  as  expected  from  the  optical  clas¬ 
sification  of  this  galaxy  as  Seyfert  1 .  The  broad  component  has  a 
line  width  ofFWHM  ~  3100  km  s_1  and  a  flux  of  5  x  10~14  ergs 
s  1  cm~2.  The  narrow  component  has  a  line  width  ofFWHM  ~ 
600  km  s_1  and  a  flux  of  1  x  10~14  ergs  s_1  cnT2. 

5.  DISCUSSION 

In  this  section  we  first  focus  on  the  investigations  of  the  en¬ 
ergy  sources  of  non-Seyfert  ULIRGs  based  on  our  /.-band  spec¬ 
tra  (§§  5. 1-5.5)  and  then  combine  these  spectra  with  data  at  other 
wavelengths  and  compare  these  sources  with  the  properties  of 
Seyfert  ULIRGs  (§§  5.6-5. 8). 

5.1.  The  Detected  Modestly  Obscured  Starbursts 

The  detection  of  the  3.3  /tm  PAH  emission  clearly  indicates 
the  presence  of  starbursts  in  the  majority  of  ULIRGs.  Since  dust 


extinction  in  the  L  band  is  about  0.06  times  as  large  as  in  the 
optical  T  band  (A  =  0.6  //m;  Rieke  &  Lebofsky  1985;  Lutz  et  al. 
1996),  the  attenuation  of  3.3  /im  PAH  emission  with  dust  extinc¬ 
tion  of  A  v  ~  15  mag  is  less  than  1  mag.  Thus,  the  observed  3 .3  /tm 
PAH  emission  luminosities  can  be  used  to  determine  the  abso¬ 
lute  magnitudes  of  modestly  obscured  (Ay  <  15  mag)  starburst 
activity. 

At  Arest  ~  3.3  /mi,  Pfr>  emission  is  present,  superposed  on  the 
3.3  /im  PAH  emission.  The  relative  contribution  from  this  Pfb 
emission  line  is  expected  to  be  the  highest  in  Seyfert  1  ULIRGs 
because  (1)  broad  components  from  the  AGNs  are  unattenuated 
and  (2)  3.3  /tm  PAH  emission  is  weak  compared  to  starburst- 
dominated  ULIRGs.  For  the  Seyfert  1  ULIRGs  Mrk  1014  and 
IRAS  07599+6508,  we  can  estimate  the  contamination  by  the  Pfr> 
emission  using  the  measured  Pa  a  emission  line  fluxes  (this  pa¬ 
per;  Taniguchi  et  al.  1994).  For  high-transition  broad  and  narrow 
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emission  lines,  such  as  Paa  and  Pf<5,  case  B  is  applicable 
( Rhee  &  Larkin  2000),  and  so  we  adopt  a  Pf<5-to-Paa  flux  ratio  of 
0.023  (Wynn- Williams  1984).  The  estimated  Pfr)  fluxes,  includ¬ 
ing  both  broad  and  narrow  components,  are  more  than  a  factor  of 
20  smaller  than  the  measured  3.3  /tm  PAH  flux  (Mrk  1014)  or  its 
upper  limit  (IRAS  07599+6508)  in  Table  3.  For  another  Seyfert 
1  UL1RG  Mrk  23 1 ,  we  find  that  the  Pfr  flux,  estimated  based  on 
the  Bry  (Arest  =  2.17  /tm)  flux  (Goldader  et  al.  1 995)  and  case  B 
(Pf<5/Br7  =  0.26;  Wynn- Williams  1984),  is  a  factor  of  >25 
smaller  than  the  observed  3.3  /tm  PAH  flux  ( Table  3).  For  the  re¬ 
maining  Seyfert  1  ULIRGs,  we  adopt  a  median  quasar  spectrum 
at  Arest  =  1 .8-2.0  /.tm  (Murphy  et  al.  1999)  and  an  average  spec¬ 
tral  energy  distribution  of  quasars  at  X  >  1  /tm  ( Fv  ex 
Neugebauer  et  al.  1987).  Assuming  case  B,  we  find  that  the 
equivalent  width  of  the  Pf<5  line  is  ^0.3  nm  for  a  typical  Seyfert  1 
source.  For  Seyfert  2  and  non-Seyfert  ULIRGs,  the  Pfr  con¬ 
taminations  are  expected  to  be  much  smaller  than  Seyfert  1 
ULIRGs.  Whenever  Paa  fluxes  are  available  in  the  literature 
( Veilleux  et  al.  1997,  1999b),  we  confirm  that  the  expected  Pff> 
fluxes  never  exceed  10%  of  the  observed  3.3  /tm  PAH  fluxes  or 
their  upper  limits  in  Seyfert  2  and  non-Seyfert  ULIRGs.  Thus, 
we  can  safely  assume  that  Pf<5  contamination  is  insignificant  for 
the  excess  component  at  Arest  ~  3.3  /im,  which  must  be  mostly 
ascribed  to  3.3  /tm  PAH  emission  (§  4.1). 

The  observed  3.3  /tm  PAH-to-infrared  luminosity  ratios 
(X3  jpah/A-ir)  are  summarized  in  column  (4)  of  Table  3.  The  ratios 
are  2  to  more  than  1 0  times  smaller  than  those  found  in  lower 
luminosity,  less  obscured  starbursts  (~10~3;  Mouri  et  al.  1990; 
Imanishi  2002).  Since  emission  from  the  ULIRGs’  cores  (<500  pc) 
is  covered  in  our  slit  spectra,  except  for  the  nearest  UL1RG  Arp 
220  (z  =  0.018),  the  small  A3.3pah/Air  ratios  suggest  that  the  de¬ 
tected  nuclear  starbursts  can  account  for  only  a  small  fraction  of 


the  infrared  luminosities  of  the  observed  ULIRGs.  Fischer  (2000) 
has  previously  found  that  the  6.2  /tm  PAH-to-infrared  luminos¬ 
ity  ratios  in  ULIRGs  are  smaller  by  a  similar  factor  than  in  less 
obscured  starburst  galaxies.  The  6.2  /tm  PAH  luminosities  were 
measured  with  the  large  apertures  of  the  Infrared  Space  Observa¬ 
tory  (ISO),  and  yet  they  show  a  similar  PAH  deficit:  this  indicates 
that  PAH  emission  from  the  host  galaxies,  outside  our  slit  spectra, 
is  not  energetically  important  in  the  infrared  luminosities  of  ULIRGs 
(Soifer  et  al.  2000).  Both  the  detected  3.3  /on  emission  and 

6.2  /im  PAH  emission  probe  modestly  obscured  (Av  <15  mag  or 
so)  starbursts,  presumably  at  the  outer  regions  of  the  ULIRGs’ 
cores.  The  energy  sources  responsible  for  the  bulk  of  the  infrared 
emission  of  these  ULIRGs  must  be  deeply  buried  (A  v  »  1 5  mag). 

5.2.  Buried  AGNs  with  Very  Weak  Starbursts 

The  energetically  dominant,  deeply  buried  energy  sources  at 
the  cores  of  non-Seyfert  ULIRGs  can  be  either  buried  AGNs  or 
very  compact  starbursts  or  both.  Since  buried  AGNs  can  produce 
large  infrared  dust  emission  luminosities,  without  PAH  emission 
(see  §  1),  the  small  observed  A3.3pah/Air  values  are  naturally 
explained.  Very  highly  obscured  starbursts  can  also  explain  the 
small  observed  A3.3pah/Air  ratios  because  flux  attenuation  of  the 

3.3  /tm  PAH  emission  may  be  significant,  while  that  of  longer 
wavelength  infrared  emission  (8-1000  /tm)  may  not.  These  two 
scenarios  are  difficult  to  differentiate  based  on  the  absolute 
luminosities  of  PAH  emission  but  can  be  distinguished  by  the 
equivalent  width  of  emission  or  absorption  features.  In  a  pure 
normal  starburst,  where  H  n  regions,  molecular  gas,  and  photo¬ 
dissociation  regions  are  spatially  mixed,  the  equivalent  width  of 
the  3.3  f.i in  PAH  emission  feature  is  less  affected  by  dust  ex¬ 
tinction  (§  1),  and  so  the  observed  values  should  be  close  to  the 
intrinsic  one  (EW3 3PAh  ~  100  nm;  Moorwood  1986;  Imanishi 
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Fig.  2. — Infrared  Z-band  (20bS  =  2.8— 4.1  /tm)  spectra  of  ULIRGs  classified  optically  as  H  n  regions.  Units  and  symbols  are  the  same  as  Fig.  1. 


&  Dudley  2000),  regardless  of  dust  extinction  of  the  starburst. 
On  the  other  hand,  an  AGN  produces  a  PAH-fr ee  featureless  con¬ 
tinuum,  with  strong  absorption  features  when  the  AGN  is  highly 
dust  obscured.  If  this  emission  contributes  significantly  to  an 
observed  Z-band  flux,  then  the  EW3.3PAH  value  should  be  smaller 
than  in  normal  starburst  galaxies.  The  EW3.3PAH  values  in  star- 
bursts  have  an  average  value  of  EW3  jpah  ~  100  nm,  with  some 
scatter,  but  never  become  lower  than  40  nm  (Moorwood  1986). 
Thus,  if  we  adopt  EW3.3PAH  ~  40  nm  as  strong  signatures  of  a 
significant  contribution  from  the  AGN ’s  PAH-free  continuum,  then 
we  can  argue  that  the  three  LINER  ULIRGs  IRAS  08572+391 5N  W, 
IRAS  12127— 1412NE,  and  IRAS  17044+6720  possess  power¬ 
ful  AGNs.  All  of  them  display  strong  dust  absorption  features  at 
3.1  /im  and/or  3 .4  gm,  suggesting  that  the  AGNs  are  deeply  buried. 

Kim  et  al.  (2002)  estimated  the  nuclear  K'-band  magnitudes 
within  the  central  4  kpc  diameter  (K'4)  for  the  ULIRGs  in  the  IRAS 


1  Jy  sample  (their  Table  3).  We  combine  these  data  with  the  Z-band 
magnitudes  derived  from  our  slit  spectra  (Zspec)  and  find  that  the 
K\  -  Zspec  colors  of  the  three  ULIRGs  IRAS  08572+3915NW, 
IRAS  12127-1412NE,  and  IRAS  17044+6720  occupy  the  larg¬ 
est  (reddest)  three  values.  Hot  dust  emission  heated  by  an  AGN 
produces  redder  K  —  L  colors  than  normal  stellar  emission  and 
so  the  K  —  Z  colors  of  AGNs  are  usually  redder  than  starbursts 
(Willner  et  al.  1984;  Alonso-Herrero  et  al.  2003;  Imanishi  & 
Alonso-Herrero  2004;  Imanishi  &  Wada  2004).  The  red  K\  — 
ZSpec  colors  again  support  the  scenario  that  powerful  buried 
AGNs  are  present  and  produce  strong  hot  dust  emission  in  these 
ULIRGs. 

5.3.  Buried  AGNs  with  Coexisting  Strong  Starbursts 

Based  on  the  EW3.3PAH  values,  the  presence  of  buried  AGNs  is 
strongly  indicated  in  the  three  LINER  ULIRGs  discussed  above. 
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Fig.  2. — Continued 


This  was  possible  primarily  because  the  3.3  /./,m  PAH  emission 
expected  from  starbursts  is  exceptionally  weak  in  these  ULIRGs. 
Hence,  these  ULIRGs  may  be  classified  as  almost  pure  buried 
AGNs.  However,  the  nuclei  of  ULIRGs  are  generally  extremely 
rich  in  molecular  gas  (Sanders  &  Mirabel  1996),  which  not  only 
fuels  the  existing  supermassive  black  holes  and  thereby  in¬ 
creases  the  luminosities  of  AGNs  but  is  also  likely  to  fuel  star 
formation.  Hence,  it  is  expected  that,  if  powerful  AGNs  are  pres¬ 
ent  in  the  majority  of  non-Seyfert  ULIRGs,  they  will  coexist  with 
starbursts.  It  is  therefore  important  to  look  for  the  signatures 
of  powerful  buried  AGNs  in  the  cores  of  starburst-containing 
ULIRGs.  To  detect  the  AGN  signatures  in  these  composite  cores 
requires  more  careful  analysis  of  the  spectra  than  for  the  almost 
pure  buried  AGNs  (Imanishi  et  al.  2001;  Soifer  et  al.  2002; 
Imanishi  &  Maloney  2003;  Spoon  et  al.  2004).  Even  if  the  in¬ 
trinsic  luminosities  of  a  buried  AGN  and  surrounding  less  ob¬ 


scured  starbursts  are  similar,  the  flux  from  the  buried  AGN  will 
be  more  highly  attenuated  by  dust  extinction.  In  addition  to  this 
extinction,  when  the  broad  3.1  /im  H20  ice  absorption  feature  is 
present  (Spoon  et  al.  2000;  Imanishi  &  Maloney  2003),  the  AGN 
flux  is  suppressed  even  more  severely.  In  this  case,  the  3.3  ^m 
PAH  emission  is  not  diluted  significantly  by  the  PAH-free  con¬ 
tinuum  produced  by  the  AGN,  and  so  the  small  EW^jpah  diag¬ 
nostic  can  no  longer  be  used. 

To  determine  whether  the  remaining  ULIRGs’  cores  consist  of 

(1)  modestly  obscured  starbursts  and  deeply  buried  starbursts  or 

(2)  modestly  obscured  starbursts  and  deeply  buried  AGNs,  we 
use  the  absolute  optical  depth  values  of  dust  absorption  features 
found  in  the  L-band  spectra.  As  mentioned  in  §  1 ,  these  values 
can  be  used  to  distinguish  whether  the  energy  sources  are  spatially 
well  mixed  with  dust  or  are  more  centrally  concentrated  than  the 
dust.  Here  we  discuss  this  point  more  quantitatively. 
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In  a  mixed  dust/source  geometry,  flux  attenuation  as  a  func¬ 
tion  of  optical  depth  is 

1  -  e~Tl 

i(n)  =  i0 - ,  (i) 

n 


where  70  is  the  unattenuated  intrinsic  flux  and  is  the  optical 
depth  at  each  wavelength;  this  takes  different  values  inside  and 
outside  the  absorption  features.  The  3.1  /im  H20  ice  absorption 
feature  is  observed  if  dust  grains  are  located  deep  inside  molecular 
gas,  are  shielded  from  ambient  U  V  radiation,  and  are  covered  with 
an  ice  mantle  ( Whittet  et  al.  1988).  The  feature  is  not  observed  in 
the  ditfuse  interstellar  medium  outside  molecular  clouds,  since  ice 
mantles  cannot  survive  exposure  to  the  ambient  UV  radiation. 
Thus,  its  optical  depth  reflects  the  total  column  density  of  ice- 
covered  dust  grains  in  molecular  gas,  in  front  of  the  continuum- 
emitting  energy  sources.  The  observed  optical  depth  of  the  3 . 1  /mi 
H20  ice  absorption  feature  is,  by  definition,  dependent  on  the  dif¬ 
ference  of  flux  attenuation  between  the  ice  absorption  feature  and 
the  3-4  /im  continuum  outside  the  feature.  Imanishi  &  Maloney 
(2003)  have  shown  that,  in  a  mixed  dust/source  geometry, 
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where/is  the  fraction  of  dust  that  is  covered  with  an  ice  mantle, 
and  it  is  assumed  that  the  T31 /Horatio  (the  abundance  of  H20  ice 
relative  to  dust  grains  in  the  regions  where  an  ice  mantle  can 


survive)  is  similar  to  that  in  Galactic  molecular  clouds  (Smith 
et  al.  1993;  Tanaka  et  al.  1990;  Murakawa  et  al.  2000). 

Stars  are  formed  within  the  high-density  cores  of  molecular 
clouds.  If  we  observe  such  molecular  clouds  individually,  then 
the  energy  sources  should  be  more  centrally  concentrated  than 
the  gas  and  accompanying  dust.  However,  starburst  galaxies 
consist  of  many  such  stars  and  the  energy  sources  will  be  distri¬ 
buted  throughout  the  volume  of  the  starburst  (Soifer  et  al.  2000; 
Siebenmorgen  et  al.  2004).  Hence,  for  observations  of  entire 
starburst  nuclei,  the  stellar  energy  sources  and  dust  will  be  rel¬ 
atively  well  mixed  spatially  (Fig.  la).  In  fact,  detailed  studies  of 
the  prototypical  starburst  galaxy  M82  strongly  suggest  that  the 
stars  of  the  burst  and  the  obscuring  dust  are  spatially  well  mixed 
(Puxley  1991;  McLeod  et  al.  1993;  Forster  Schreiber  et  al.  2001). 
Thus,  it  seems  quite  reasonable  to  regard  that  starbursts  can  be 
represented  by  this  mixed  dust/source  geometry  model. 

The  observed  T3.1  depends  on+Fand  the  Rvalue,  the  fraction 
of  dust  that  is  covered  with  an  ice  mantle.  If  we  assume  that  the 
Rvalues  at  the  cores  of  ULIRGs  are  similar  to  that  in  the  well- 
studied  prototype  starburst  galaxy  M82  (/  ~  30%;  Imanishi  & 
Maloney  2003),  then  r3  1  will  be  <0.3  for  any  Av  value  in  a 
mixed  dust/source  geometry  (Imanishi  &  Maloney  2003).  In 
reality,  the  cores  of  ULIRGs  are  estimated  to  have  much  higher 
surface  brightnesses  than  M82  (Soifer  et  al.  2000),  implying  a 
much  larger  ambient  radiation  field.  In  a  mixed  dust/source 
geometry,  the  radiation  environment  in  the  cores  of  ULIRGs  will 
therefore  be  much  harsher  and  it  will  be  much  more  difficult  for 
ice  mantles  to  survive  than  in  M82.  In  this  geometry,  the  Rvalues 
in  ULIRGs’  cores  must  be  smaller  than  the  assumed  value  of 
/  =  30%,  and  thus  the  upper  limit  of  r3.i  =  0.3  should  be  robust. 
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Fig.  3. — Infrared  L-band  (20bs  =  2.8— 4.1  fim)  spectra  of  optically  unclassified  ULIRGs.  Units  and  symbols  are  the  same  as  Fig.  1. 
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Fig.  4. — Infrared  L-band  (/l0bS  =  2.8— 4.1  /xm)  spectra  of  ULIRGs  classified  optically  as  Seyfert  2  galaxies.  Units  and  symbols  are  the  same  as  Fig.  1. 

The  r 3  i  value  can  become  arbitrarily  large  with  increasing  Ay 
provided  that /  is  nonzero.  In  this  centrally  concentrated  energy- 
source  geometry,  once  the  shielding  column  density  of  dust  and 
molecular  gas  around  the  central  energy  source  becomes  suffi¬ 
ciently  large,  dust  grains  at  a  given  distance  from  the  center  can 
be  covered  with  ice  mantles  (/  >  0),  and  so  a  large  T3.1  value  can 
be  produced.  Therefore,  observing  larger  7-3.1  values  than  the 
limit  for  a  mixed  dust/source  geometry  (0.3)  is  another  inde¬ 
pendent  means  to  determine  the  geometry  of  ULIRGs’  cores  and 
therefore  the  nature  of  the  central  energy  sources. 

Aside  from  the  three  LINER  ULIRGs  with  strong  buried  AGN 
signatures  discussed  above  (IRAS  08572+39 15NW,  IRAS 
12127— 1412NE,  and  IRAS  17044+6720),  14  LINER,  three  Hu, 
and  two  optically  unclassified  ULIRGs’  nuclei  have  observed 


On  the  other  hand,  m  a  buried  AGN,  the  energy  source  is  very 
compact  and  is  more  centrally  concentrated  than  the  surrounding 
dust  (Fig.  lb,  left  panel).  In  this  geometry,  a  simple  foreground 
screen  dust  model  is  applicable  and  the  flux  attenuation  is 

I(n)  =  ke~TK  (3) 

The  T3.i  value  is  expressed  as 

q  7"cont 
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=  0.0  6Avf.  (4) 
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t3  i  values  larger  than  the  threshold  ( Table  4).  Note  that  these  ob¬ 
served  values  are  contaminated  by  weakly  obscured  starburst  emis¬ 
sion.  Hence,  these  observed  optical  depths  provide  only  lower 
limits  to  the  true  r3  4  values  for  the  central  buried  energy  sources. 
For  these  ULIRGs’  nuclei,  a  centrally  concentrated  energy-source 
geometry  provides  a  natural  explanation  for  the  large  observed 
t3  i  values. 

There  is  one  case  in  which  a  mixed  dust/source  geometry  for 
the  cores  can  exceed  the  threshold  value  of  0.3.  If  the  host  gal¬ 
axies  are  viewed  close  to  edge-on  and  a  large  column  density  of 
ice-covered  dust  is  present  in  the  interstellar  medium  of  the  host 
galaxies,  exterior  to  the  starburst  core,  then  a  foreground  screen 
model  is  applicable  and  r31  can  be  large,  just  as  for  a  centrally 
concentrated  energy  source.  However,  the  presence  of  ice-covered 
dust  requires  that  the  dust  be  located  deep  inside  molecular  gas 
( Whittet  et  al.  1988).  Strong  3.1  /mi  H20  ice  absorption  is  found 
in  the  dense  cores  of  molecular  gas  (Murakawa  et  al.  2000),  and 
such  dense  molecular  gas,  probed  by,  e.g.,  HCN  (./  =  1-0) 
emission  at  v  —  88.632  GHz  (1  =  3.3848  mm)  in  the  millimeter 
wavelength  range,  is  normally  seen  only  in  the  nuclei  of  ULIRGs 
(Imanishi  et  al.  2005).  It  is  thus  very  likely  that  the  strong  ob¬ 
served  H20  ice  absorption  largely  arises  in  the  cores  of  the 
ULIRGs  rather  than  in  the  extended  interstellar  medium  of  the  host 
galaxies.  The  large  observed  H20  ice  absorption  optical  depths  can 
thus  be  taken  as  some  evidence  of  a  centrally  concentrated  energy- 
source  geometry  at  the  cores  of  these  ULIRGs. 

Soifer  et  al.  (2000)  studied  the  cores  of  seven  nearby  ULIRGs 
and  estimated  that,  in  the  majority  (six  out  of  seven)  of  these 
ULIRGs’  cores,  the  emission  surface  brightnesses  are  very  high 
(~1OI3Z0  kpc  2)  and  close  to  the  maximum  observed  values  seen 
in  the  cores  of  Galactic  H  n  regions,  even  if  the  energy  sources  are 
uniformly  distributed  over  the  cores  (<500  pc  in  size).  Provided 
that  the  emission  properties  of  the  six  ULIRGs’  cores  (Soifer  et  al. 
2000)  are  representative  of  the  majority  of  nearby  ULIRGs,  sur¬ 
face  brightnesses  of  the  centrally  concentrated  energy  sources 
would  be  even  higher,  making  exceptionally  centrally  concen¬ 
trated  starbursts  less  likely,  if  not  completely  ruled  out.  AGNs  are 
more  natural  explanations  for  the  centrally  concentrated  high  sur¬ 
face  brightness  energy  sources  at  the  cores  of  ULIRGs.  We  thus 
take  the  centrally  concentrated  energy  sources  as  tentative  evi¬ 
dence  for  buried  AGNs. 

We  can  apply  the  same  argument  concerning  the  geometry  of 
the  energy  sources  and  dust  to  the  absolute  optical  depth  of  the 
3.4  pm  carbonaceous  dust  absorption  feature  (r3  4).  This  ab¬ 
sorption  feature  is  found  if  a  source  is  obscured  by  bare  dust 
grains  (Pendleton  et  al.  1994;  Imanishi  et  al.  1996;  Rawlings  et  al. 
2003),  but  absent  if  the  absorbing  dust  is  ice  covered  ( Mennella 


et  al.  2001).  Thus,  the  optical  depth  of  the  3.4  pm  absorption  fea¬ 
ture  traces  the  column  density  of  bare  carbonaceous  dust  grains. 
Imanishi  &  Maloney  (2003)  have  quantitatively  demonstrated 
that  r3  4  cannot  be  larger  than  0.2  in  a  mixed  dust/source  geom¬ 
etry.  Two  LINER  ULIRGs,  IRAS  08572+3915NW  and  IRAS 
12127-  1412NE,  show  r3.4  >  0.2  (Table  4).  These  ULIRGs  have 
already  been  classified  as  buried  AGNs  based  on  their  EW3  3PAH 
values  (§  5.2).  For  them,  additional  support  for  powerful  buried 
AGNs  is  provided  here  from  the  geometry  implied  by  the  depth  of 
their  3 .4  pm  absorption  features. 

For  the  remaining  ULIRGs  where  (1)  the  EW3 ,3pah  values  are 
in  the  possible  scattered  range  of  starburst  galaxies  (>40  nm)  and 
(2)  no  strong  dust  absorption  features  are  observed,  there  are 
no  obvious  signatures  of  buried  AGNs.  Their  spectra  can  be  ex¬ 
plained  by  either  of  the  following  scenarios:  (1)  the  cores  are 
dominated  by  starbursts  with  a  mixed  dust/source  geometry,  or 
(2)  the  cores  do  possess  powerful  buried  AGNs,  but  the  emission 
from  the  AGNs  is  so  highly  attenuated  that  its  contribution  to  the 
observed  L-band  fluxes  is  not  significant.  We  have  no  way  to  dis¬ 
tinguish  between  these  two  scenarios,  and  their  energy  sources 
remain  undetermined  by  our  /.-band  spectroscopic  method.  There 
are  some  known  examples  (e.g.,  NGC  4418)  that  show  no  obvious 
buried  AGN  signatures  in  the  L-band  spectra  (Imanishi  et  al.  2004) 
but  almost  certainly  contain  powerful  AGNs  based  on  longer 
wavelength  data  (Dudley  &  Wynn-Williams  1997;  Spoon  et  al. 
2001).  Observations  at  wavelengths  with  lower  dust  extinction  are 
necessary  to  unveil  possible  buried  AGNs  that  are  undetectable 
even  in  the  L  band. 

5.4.  Dust  Extinction  and  Dereddened  AGN  Luminosities 

In  a  centrally  concentrated  energy-source  geometry,  dust  at 
a  temperature  of  1000  K,  which  is  close  to  the  innermost  dust 
sublimation  radius,  produces  emission  with  a  peak  at  X  ~  3  pm 
( Fig.  lb,  right  panel),  assuming  approximately  blackbody  emis¬ 
sion.  Dust  grains  at  larger  radii,  with  lower  temperatures,  show 
emission  peaks  at  longer  infrared  wavelengths.  The  dust  ex¬ 
tinction  toward  the  3  pm  continuum-emitting  region,  Av( 3  pm), 
can  be  estimated  from  the  r3  4  (bare  dust)  and  r34  (ice-covered 
dust)  values  and  can  provide  a  better  estimate  for  dust  extinction 
toward  the  buried  AGN  itself  than  the  estimate  using  longer 
infrared  wavelengths,  which  only  probes  dust  extinction  toward 
the  outer,  cooler  emission  regions  (Fig.  lb,  right  panel).  If  we 
assume  that  the  optical  depth  of  these  dust  absorption  features, 
relative  to  Ay,  is  similar  to  the  Galactic  diffuse  interstellar  me¬ 
dium  (t3.4/Hk  =  0.004-0.007,  Pendleton  et  al.  1994;  r^.i/Ay  = 
0.06,  Smith  et  al.  1993;  Tanaka  et  al.  1990;  Murakawa  et  al. 
2000),  and  if  we  adopt  a  foreground  screen  model  that  should  be 
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Fig.  5. — Infrared  L-band  (20bS  =  2.8— 4.1  fim)  spectra  of  ULIRGs  classified  optically  as  Seyfert  1  galaxies.  Units  and  symbols  are  the  same  as  Fig.  1. 


applicable  to  the  centrally  concentrated  energy-source  geometry, 
then  the  estimated  Av  values  are  summarized  in  column  (4)  of 
Table  4. 

For  the  three  ULIRGs  with  very  small  EW3  3PAH  values,  IRAS 
08572+39 15NW,  IRAS  12127-1412NE,  and  IRAS  17044+ 
6720,  the  /.-band  emission  is  concluded  to  predominantly  come 
from  hot  dust  emission  heated  by  the  buried  AGN.  We  can  esti¬ 
mate  the  dereddened  luminosity  of  the  AGN -heated  dust  emis¬ 
sion  fairly  straightforwardly  for  these  sources.  Assuming  that 
dust  extinction  atZ,  relative  to  V,  is  similar  to  the  Galactic  diffuse 
interstellar  medium  ( AJAy  ~  0.06;  Rieke  &  Lebofsky  1985; 
Lutz  et  al.  1996),  the  dereddened  AGN  luminosities  at  Arest  = 
3-4  /im  are  >50  x  1045,  >10  x  1045,  and  (0.7-2)  x  1045  ergs  s-1, 
respectively,  for  IRAS  08572+3915NW,  IRAS  12127— 


1412NE,  and  IRAS  17044+6720.  The  dereddened  AGN  lumi¬ 
nosities  could  quantitatively  account  for  the  bulk  of  the  infrared 
luminosities  of  these  ULIRGs  (Lir  ~  5  x  1045  ergs  s_1).  For 
IRAS  08572+3915NW  and  IRAS  12127-1412NE,  the  dered¬ 
dened  AGN  luminosities  significantly  exceed  the  observed  in¬ 
frared  luminosities.  A  similar  excess  was  found  in  the  ULIRG 
Superantennae  (Risaliti  et  al.  2003),  and  the  excess  implies  that 
the  Tt,j/Al  ratios  in  some  ULIRGs’  cores  are  larger  than  the 
ratios  found  in  the  Galactic  interstellar  medium. 

For  the  remaining  starburst/buried-AGN  composite  ULIRGs 
(large  EW3  3PAH  values  with  >40  nm  and  strong  3.1  /im  H20  ice 
absorption  with  r3.i  >  0.3),  no  detectable  3.4  /im  absorption 
features  are  found.  The  optical  depths  of  the  3.1  and  3.4  /im 
absorption  features  probe  the  column  density  of  ice-covered  and 
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bare  dust  grains,  respectively,  but  it  seems  strange  that  there  are 
no  bare  dust  grains  in  the  nuclear  regions  of  these  ULIRGs.  The 
3.3  //ni  PAH  emission  feature  is  often  accompanied  by  a  subpeak 
at  3.4  /im  (Tokunaga  et  al.  1991;  Imanishi  &  Dudley  2000),  and 
this  subpeak  may  dilute  the  3.4  /im  dust  absorption  feature  at 
the  same  wavelength.  In  fact,  some  ULIRGs’  spectra  with  strong 
3.3  pm  PAH  emission  show  subpeaks  at  Arest  =  3.4  /im  in  Fig¬ 
ures  1,2,  and  3  (e.g.,  IRAS  09039+0503,  IRAS  13509+0442,  and 
IRAS  14197+0813).Forthose ULIRGs withEW3.3PAH  >  40 nm, 
the  estimated  Ay  values  (Table  4,  col.  [4])  should  be  taken  as  only 
lower  limits  for  dust  extinction  toward  the  inner  3  /tm  continuum- 
emitting  regions  at  the  cores.  Although  quantitative  values  are 
highly  uncertain  for  these  sources,  the  cores  are  usually  assumed 
to  dominate  the  infrared  luminosities  of  nearby  ULIRGs  (Soifer 
et  al.  2000),  and  the  cores  of  these  composite  ULIRGs  are  sug¬ 
gested  to  have  centrally  concentrated  energy  sources,  most  likely 
buried  AGNs.  Thus,  the  buried  AGNs  in  these  ULIRGs  could 
account  for  the  bulk  of  their  observed  infrared  luminosities. 

5.5.  Dependence  on  Optical  Spectral  Type 

In  total,  /.-band  spectra  of  27  LINER  and  13  H  n  ULIRGs’ 
nuclei  have  been  obtained.  Three  LINER  ULIRGs’  and  no  H  n 
ULIRGs’  nuclei  show  clear  evidence  for  almost  pure  buried 
AGNs  (§  5.2).  Additionally,  14  LINER  and  three  H  n  ULIRGs’ 
nuclei  have  tentative  evidence  for  powerful  buried  AGNs,  in 
addition  to  modestly  obscured  starbursts  (§  5.3).  Combining 
these  values,  17  out  of  27  LINER  ULIRGs’  nuclei  (63%)  and  3 
out  of  13  H  ii  ULIRGs  (23%)  display  some  degree  of  evidence 
for  powerful  buried  AGNs.  The  observed  LINER  ULIRGs 
comprise  a  complete  sample  (§  2).  The  observed  H  n  ULIRGs 
also  constitute  a  complete  sample  in  the  right  ascension  range 
10-22  hr  (§  2).  The  fraction  of  sources  with  signatures  of  buried 
AGNs  is  significantly  higher  in  LINER  ULIRGs  than  H  n 
ULIRGs.  This  suggests  that  a  larger  fraction  of  LINER  ULIRGs 
may  have  powerful  buried  AGNs  than  H  n  ULIRGs. 

Some  ULIRGs  have  spectra  of  poor  signal-to-noise  ratios,  and 
the  presence  of  strong  dust  absorption  features  at  3. 1  pm  (r3.i  > 
0.3)  or  3.4  /tm  (r3.4  >  0.2)  may  be  missed.  Strong  3.1  pm  H20 
ice  absorption  may  be  present,  but  be  missed  due  to  poor  signal- 
to-noise  ratios,  in  four  LINER  ULIRGs’  nuclei  (IRAS  12112+ 
0305NE,  IRAS  12112+0305SW,  IRAS  14252-1550W,  and 
IRAS  16468+5200E)  in  Figure  1  and  one  H  n  ULIRG’s  nucleus 
(IRAS  21208-05 19N)  in  Figure  2.  Two  LINER  ULIRGs 
(IRAS  10485-1447  and  IRAS  16090-0139)  may  have  large 
r3.4  values.  Even  if  only  the  H  n  UL1RG  IRAS  21208-0519N  is 
included  as  sources  with  tentative  buried  AGN  signatures,  the 


AGN  fraction  in  H  n  ULIRGs’  nuclei  is  only  31%  (4/13),  sig¬ 
nificantly  below  the  fraction  for  LINER  ULIRGs’  nuclei. 

In  a  pure  buried  AGN,  XDRs,  dominated  by  low-ionization 
species  (Maloney  et  al.  1996),  should  develop  and  aLINER-type 
optical  spectrum  is  expected  (Imanishi  et  al.  2001),  rather  than  a 
typical  Seyfert  spectrum  originating  in  the  NLRs.  The  discovery 
of  pure  powerful  buried  AGNs  with  no  detectable  starbursts  in 
the  two  LINER  ULIRGs  (IRAS  08572+39 15NW  and  IRAS 
12127— 1412NE)  may  be  explained  by  this  scenario. 

For  the  remaining  LINER  ULIRGs  with  coexisting  modestly 
obscured  starbursts,  the  LINER  optical  spectra  can  be  due  to 
XDRs  and/or  shocks  driven  by  the  starbursts  (Taniguchi  et  al. 
1999;  Lutz  et  al.  1999).  In  the  case  of  XDRs,  even  though  the 
central  AGN  itself  is  deeply  buried,  emission  from  the  surround¬ 
ing  XDRs  can  be  relatively  unobscured.  Although  the  origin  of 
the  LINER  optical  spectra  in  these  ULIRGs  is  unclear,  the  higher 
detection  rate  of  powerful  buried  AGNs  in  LINER  ULIRGs 
suggests  that  emission  from  XDRs  may  be  important  in  the 
optical. 

5.6.  Comparison  with  VLA  Radio  Observations 

Nagar  et  al.  (2003)  performed  VLA  observations  at  v  —  15  GHz 
(A  =  2  cm),  with  150  mas  spatial  resolution,  of  a  large  sample  of 
ULIRGs  in  the  same  IRAS  1  Jy  sample  and  searched  for  compact 
radio  core  emission,  another  good  indicator  of  an  AGN.  They 
concluded  that  the  detection  rate  of  the  radio  cores  in  LINER 
ULIRGs  (~74%)  was  substantially  higher  than  in  H  n  ULIRGs 
(^14%).  Our  higher  detection  rate  of  buried  AGN  signatures  in 
LINER  ULIRGs  than  in  H  n  ULIRGs  is  both  qualitatively  and 
quantitatively  similar  to  their  conclusion. 

Table  3  (cols.  [6]  and  [7])  summarizes  information  on  which 
individual  ULIRGs  have  detectable  (or  nondetectable)  radio 
cores  (Nagar  et  al.  2003),  in  comparison  with  the  buried  AGN 
signatures  from  our  Z-band  spectroscopic  method.  Among  the 
three  ULIRGs  with  clear  buried  AGN  signatures,  two  sources, 
IRAS  08572+3915NW  and  IRAS  17044+6720,  were  observed 
with  VLA  and  both  (100%  =  2/2)  show  core  radio  emission. 
For  ULIRGs  with  tentative  evidence  for  buried  AGNs,  eight 
LINER,  three  H  n,  and  one  optically  unclassified  ULIRGs’  nuclei 
have  VLA  data,  and  five  LINER,  two  H  n,  and  one  unclassified 
ULIRGs’  nuclei  have  detectable  VLA  radio  cores.  Hence,  the 
detection  rate  is  67%  [=(5  +  2  +  l)/(8  +  3  +  1)].  For  ULIRGs 
without  buried  AGN  signatures  in  the  Z-band  spectra,  VLA  cores 
were  detected  in  8  out  of  17  sources  (47%).  The  higher  detection 
rate  of  the  VLA  cores  in  ULIRGs  that  show  a  higher  level  of 
buried  AGN  evidence  in  the  Z  band  suggests  that  our  Z-band 


TABLE  3 

Properties  of  the  Nuclear  3.3  fim  PAH  Emission  and  Buried  AGN  Signatures 


Object 

(i) 

(xicr 

./Lpah 

14  —1 

ergs  s 

(2) 

cm  2) 

^3.3PAH 

(xl041ergs  s-1) 

(3) 

^3.3PAH^IR 

(xlO-3) 

(4) 

Rest  EW3.3PAh 
(nm) 

(5) 

AGN  Evidence 

L  Band  Radio 

(6)  (7) 

IRAS  00188-0856 . 

3.5 

13 

0.15 

50 

O 

O  (0-6) 

IRAS  03250+1606 . 

3.0 

12 

0.25 

80 

0 

X 

IRAS  08572+39 15NW . 

<3.5a 

<2.5a 

<0.05a 

<5a 

© 

O  (0.8) 

IRAS  09039+0503 . 

1.5 

5.5 

0.10 

95 

0 

X 

IRAS  091 16+0334W . 

3.0 

15 

0.25 

75 

0 

X 

IRAS  09539+0857 . 

0.9 

3.5 

0.07 

65 

X 

0(1-0) 

IRAS  10378+1108 . 

1.0 

4.5 

0.06 

40 

0 

O  (0.8) 

IRAS  10485-1447 . 

2.0 

7.0 

0.10 

80 

0 

IRAS  10494+4424 . 

3.5 

6.5 

0.10 

110 

0 

O  (0-4) 

IRAS  11095-0238 . 

2.5 

6.0 

0.10 

150 

X 

IRAS  121 12+0305NE . 

3.5 

4.0 

0.05 

100 

X 

O  (0.8) 

IRAS  121 12+0305SW . 

3.5 

4.0 

0.05 

160 

X 

X 

IRAS  12127— 1412NE . 

0b 

0b 

0b 

0b 

© 

IRAS  12359— 0725N . 

2.0 

9.5 

0.20 

75 

0 

IRAS  14252— 1550E . 

0.7 

3.5 

0.06 

60 

0 

IRAS  14252— 1550W . 

0.7 

3.5 

0.06 

60 

X 

IRAS  14348- 1447NE . 

1.5 

2.0 

0.03 

70 

0 

IRAS  14348- 1447SW . 

3.5 

5.0 

0.07 

110 

X 

IRAS  15327+2340  (Arp  220) . 

25 

1.5 

0.03 

80 

X 

O  (0-2) 

IRAS  16090-0139 . 

3.0 

13 

0.10 

75 

O 

IRAS  16468+5200E . 

1.0 

5.5 

0.10 

120 

X 

X 

IRAS  16487+5447W . 

2.5 

6.0 

0.10 

70 

O 

O  (0-7) 

IRAS  17028+5817W . 

4.0 

10 

0.20 

120 

0 

O  (0-9) 

IRAS  17044+6720 . 

2.0 

8.5 

0.15 

10 

© 

O  (1-3) 

IRAS  21329-2346 . 

1.0 

4.0 

0.08 

50 

0 

IRAS  23234+0946W . 

2.5 

10 

0.20 

75 

X 

O  (0-4) 

IRAS  23327+2913S . 

2.0 

5.0 

0.10 

45 

X 

X 

IRAS  10190+1322E . 

5.5 

7.0 

0.20 

95 

X 

O(i.o) 

IRAS  101 90+ 1322 W . 

1.0 

1.5 

0.04 

50 

X 

X 

IRAS  11387+4116 . 

1.5 

8.0 

0.15 

75 

X 

X 

IRAS  11506+1331  . 

5.5 

21 

0.25 

95 

O 

O  (0.3) 

IRAS  13443+0802NE . 

1.0 

5.0 

0.08 

75 

X 

O  (0-7) 

IRAS  13509+0442 . 

4.0 

18 

0.25 

135 

X 

X 

IRAS  13539+2920NW . 

4.5 

11 

0.25 

85 

X 

O  (0-2) 

IRAS  14060+2919 . 

7.5 

23 

0.45 

150 

X 

X 

IRAS  15206+3342 . 

4.0 

14 

0.25 

55 

X 

X 

IRAS  15225+2350 . 

1.5 

7.0 

0.15 

40 

O 

O  (1.0) 

IRAS  16474+3430S . 

5.0 

14 

0.25 

105 

O 

X 

IRAS  20414-1651 . 

2.5 

4.0 

0.05 

75 

X 

O  (0.7) 

IRAS  21208-05 19N . 

1.5 

6.0 

0.10 

100 

X 

X 

IRAS  14197+0813 . 

2.0 

7.5 

0.15 

90 

O 

Oc(?) 

IRAS  14485-2434 . 

2.5 

13.5 

0.25 

70 

O 

IRAS  05189-2524 . 

11 

4.0 

0.09 

4 

0 

IRAS  08559+1053 . 

3.5 

17 

0.25 

10 

0 

X 

IRAS  12072— 0444S . 

1.0 

3.5 

0.04 

45 

© 

IRAS  12072— 0444N . 

2.0 

8.0 

0.09 

20 

0 

IRAS  13428+5608  (Mrk  273) . 

14 

4.0 

0.09 

35 

0 

O  (0.4) 

IRAS  1345 1+1232W  (PKS  1345+12) . 

<2.5 

<7.0 

<0.10 

<5 

© 

O  (0-9) 

IRAS  15130-1958 . 

3.5 

9.0 

0.15 

10 

0 

IRAS  17179+5444 . 

1.0 

6.0 

0.09 

20 

© 

O  (0-9) 

IRAS  01572+0009  (Mrk  1014) . 

3.0 

21 

0.15 

7 

0 

O  (i.o) 

IRAS  07599+6508 . 

<7.5 

<40 

<0.35 

<4 

0 

O  (0.8) 

IRAS  11598-0112 . 

1.5 

7.0 

0.06 

4 

© 

IRAS  12265+0219  (3C  273) . 

<7.0 

<40 

<0.20 

<2 

0 

0(1-0) 

IRAS  12540+5708  (Mrk  231) . 

18 

6.0 

0.05 

2 

0 

O  (0.9) 

IRAS  15462-0450 . 

3.5 

7.0 

0.15 

10 

0 

IRAS  21219-1757 . 

1.0 

2.5 

0.05 

1 

0 

O  (1-0) 

Notes. — Col.  (1):  Object  name.  Col.  (2):  Observed  flux  of  the  3.3  fim  PAH  emission.  The  second  effective  digit,  if  smaller  than  unity,  is  shown  in  units  of  0.5. 
Col.  (3):  Observed  luminosity  of  the  3.3  fin i  PAH  emission.  Col.  (4):  Observed  3.3  fim  PAH-to-inffared  luminosity  ratio  in  units  of  10-3,  a  typical  value  for  less 
infrared-luminous  starbursts  (Mouri  et  al.  1990;  Imanishi  2002).  Col.  (5):  Rest-frame  equivalent  width  of  the  3.3  fin l  PAH  emission.  Those  for  starbursts  are  ~100  nm 
(Moorwood  1986;  Imanishi  &  Dudley  2000).  Col.  (6):  Evidence  of  an  AGN  from  our  L-band  spectroscopic  method.  ©:  strong  evidence;  0:  tentative  evidence; 
X:  no  evidence.  Col.  (7):  VLA  radio  data  by  Nagar  et  al.  (2003).  0:  core  detected;  X:  undetected;  ellipses:  no  data.  The  values  in  parentheses  are  the  core-to-total 
radio  flux  ratio  in  mJy,  taken  from  Nagar  et  al.  (2003). 

a  A  conservative  upper  limit  of  the  3.3  fim  PAH  emission  was  estimated  based  on  our  new  spectrum  and  is  larger  than  our  original  estimate  ( Imanishi  &  Dudley  2000). 
b  The  upper  limit  of  the  3.3  fim  PAH  emission  is  difficult  to  estimate  because  of  the  winding  spectral  shape.  The  observed  profiles  are  explained  by  absorption 
features,  and  so  we  adopt  the  PAH  flux  as  zero. 

c  We  classify  IRAS  14197+0813  as  a  ULIRG  with  detectable  core  emission,  based  on  col.  (9)  of  Table  1  in  Nagar  et  al.  (2003). 
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TABLE  4 

Absorption  Features  in  the  L  Band 


Object 

(1) 

Observed  r3  i 
(2) 

Observed  t3  4 
(3) 

Ay 

(mag) 

(4) 

Optical  Class 

IRAS  00188-0856 . 

1.8 

>30 

LINER 

IRAS  03250+1606 . 

0.8 

>13 

LINER 

IRAS  08572+3915NW . 

0.8 

110-200 

LINER 

IRAS  09039+0503 . 

0.8 

>13 

LINER 

IRAS  09116-03  MW . 

0.4 

>6 

LINER 

IRAS  10378+1108 . 

0.6 

>10 

LINER 

IRAS  10485+1108 . 

0.8 

>13 

LINER 

IRAS  10494+4424 . 

1.0 

>16 

LINER 

IRAS  12127— 1412NE . 

0.4 

0.35 

56-94 

LINER 

IRAS  12359— 0725N . 

0.5 

>8 

LINER 

IRAS  14252— 1550E . 

0.7 

>11 

LINER 

IRAS  14348- 1447NE . 

0.5 

>8 

LINER 

IRAS  16090-0139 . 

0.8 

>13 

LINER 

IRAS  16487+5447W . 

1.2 

>20 

LINER 

IRAS  1 7028+58 17W . 

0.6 

>10 

LINER 

IRAS  17044+6720 . 

0.15 

21-38 

LINER 

IRAS  21329-2346 . 

1.0 

>16 

LINER 

IRAS  11506+1331 . 

1.0 

>16 

H  n 

IRAS  15225+2350 . 

0.4 

>6 

H  11 

IRAS  16474+343  OS . 

0.8 

>13 

H  11 

IRAS  14197+0813 . 

0.4 

>6 

Unclassified 

IRAS  14485-2434 . 

0.4 

>6 

Unclassified 

IRAS  12072— 0444N . 

0.4 

57-100 

Sy2 

Notes. — Col.  (1):  Object  name.  Only  ULIRGs  that  show  clear  absorption  features  are  listed.  Col.  (2): 
Observed  optical  depth  of  the  3.1  /zm  H2O  ice  absorption  feature.  Since  relatively  low  plausible  continuum  levels 
are  adopted,  these  values  should  be  taken  as  conservative  lower  limits.  The  r31  values  are  derived  by  combining  a 
few  data  points  around  the  ice  absorption  peak,  and  the  statistical  errors  are  at  a  level  of  <0. 1 .  Col.  (3):  Observed 
optical  depth  of  the  3.4  /zm  bare  carbonaceous  dust  absorption  feature.  Col.  (4):  Dust  extinction  in  ^toward  the 
3-4  /zm  continuum  emission  regions,  derived  from  the  absorption  optical  depths,  assuming  the  Galactic  dust 
extinction  curve.  For  those  with  detected  3.1  /zi n  absorption,  but  with  no  measurable  r3  4  values,  the  derived  Ay 
values  should  be  taken  as  lower  limits  (see  text  in  §  5.4). 


spectroscopic  method  is  a  reliable,  if  not  perfect,  tracer  of  elusive 
buried  AGNs  in  non-Seyfert  ULIRGs. 

Our  infrared  Z-band  spectra  can  probe  dust  emission  heated  by 
an  AGN  and  thus  quantitatively  determine  the  dust  emission  lumi¬ 
nosity,  the  dominant  emission  mechanism  of  ULIRGs.  However, 
the  spectra  may  miss  the  AGN  signatures  if  the  AGN  emission  is 
highly  attenuated  by  dust  extinction.  On  the  other  hand,  the  radio 
observations  are  negligibly  affected  by  dust  extinction  and  so 
could  detect  highly  dust-obscured  AGNs  that  may  be  missed  in 
our  /.-band  spectroscopic  method.  However,  the  radio  data  trace 


synchrotron  emission  and  may  not  be  proportional  to  ULIRGs’ 
dust  emission  luminosities  (Nagar  et  al.  2003).  Thus,  the  infrared 
Z- band  and  radio  observations  are  complementary  to  each  other. 
The  fact  that  both  independent  methods  provide  a  similar  result 
suggests  that  the  higher  detection  rate  ofpowerful  buried  AGNs  in 
LINER  ULIRGs  than  H  n  ULIRGs  is  probably  real. 

5.7.  Comparison  with  Seyfert  ULIRGs 

We  compare  the  Z-band  spectral  properties  of  non-Seyfert 
ULIRGs  showing  buried  AGN  signatures  with  those  of  Seyfert 


m 


Fig.  6. — Infrared  AT-band  (20bs  =  2.07—2.5  /zm)  spectra  of  IRAS  08572+3915  andMrk  1014,  obtained  with  IRTF  SpeX.  The  abscissa  is  the  observed  wavelength  in /zm. 
The  ordinate  is  F;  in  10-15  W  m-2  /zm“ 1 .  The  solid  line  in  the  Mrk  1014  spectrum  is  the  adopted  continuum  level  to  estimate  the  Paa  emission  flux. 


134 


1MAN1SHI,  DUDLEY,  &  MALONEY 


Vol.  637 


gradient 


Av(3um) 

Av(10um) 

Av(20um) 


Fig.  7. — (a)  Geometry  of  energy  sources  and  dust  in  a  normal  starburst.  The  energy  sources  (stars)  and  dust  are  spatially  well  mixed,  (b)  Left :  Geometry  of  an  energy 
source  and  dust  in  a  buried  AGN.  The  energy  source  is  more  centrally  concentrated  than  the  surrounding  dust.  Right :  Schematic  diagram  of  a  strong  dust  temperature 
gradient  around  a  centrally  concentrated  energy  source.  Inner  dust  has  higher  temperature.  The  blackbody  radiation  from  inner  1000  K  dust,  close  to  the  dust  sublimation 
radius,  has  a  peak  at  X  ~  3  pm.  That  of  outer  300  K  dust  and  that  of  even  outer  1 50  K  dust  have  peaks  at  X  ~  10  and  ~20  pm,  respectively.  This  is  a  very  simplified  picture, 
and  in  an  actual  case,  dust  with  other  temperatures  can  also  contribute  to  the  observed  flux  at  each  wavelength  of  3, 1 0,  and  20  /im.  The  contribution  is  highly  dependent  on 
the  dust  radial  density  distribution,  but  regardless  of  any  dust  radial  density  distribution,  it  is  always  true  that  a  contribution  from  inner,  hotter  dust  emission  is  larger  at 
shorter  wavelengths  in  this  wavelength  range.  Thus,  the  relation  Ay(3  /im)  >  Ay(l0  /im)  >  Ay( 20  /im)  should  hold,  and  dust  extinction  estimated  using  3  /im  data  is  a 
better  probe  for  dust  column  density  toward  the  AGN  itself  than  longer  wavelengths.  Dust  extinction  toward  the  3  /im  continuum-emitting  region,  Ay (3  /im),  can  be 
estimated  from  T3.4  (bare  dust)  and  Ty\  (ice-covered  dust)  values.  Dust  extinction  toward  10  and  20  /im  continuum-emitting  regions,  Ay(  10  /im)  and  Ay (20  /im),  can  be 
estimated  from  the  optical  depths  of  9.7  and  18  /im  silicate  dust  absorption  features,  respectively.  These  Ay(l0  /im)  and  ,4^(20  /im)  data  will  be  available  for  all  the  non- 
Seyfert  ULIRGs  studied  in  this  paper  from  Spitzer  IRS  spectra  of  guaranteed  time  observations  and  Cycle  1  programs,  including  ours.  Comparison  of  these  values  can  be 
used  to  further  investigate  whether  the  energy  source  is  more  centrally  concentrated  than  the  surrounding  dust  (Dudley  &  Wynn-Williams  1997;  Imanishi  2000). 


ULIRGs  (i.e.,  known  AGN -possessing  ULIRGs).  The  apparent 
main  difference  between  them  is  that  the  ionizing  radiation  from 
the  putative  AGNs  in  non-Seyfert  ULIRGs  is  obscured  by  the 
surrounding  dust  along  virtually  all  lines  of  sight,  while  dust 
around  the  AGNs  in  Seyfert  ULIRGs  is  distributed  in  a  “torus,” 
and  ionizing  radiation  from  the  AGNs  can  escape  along  the  torus 
axis,  allowing  the  Seyfert  signatures  to  be  detectable  in  the  op¬ 
tical.  In  Table  3,  all  Seyfert  ULIRGs’  nuclei  show  EW3.3PAH 
smaller  than  ~40  nm.  The  EW3  3PAh  values  of  Seyfert  1  ULIRGs 
are  systematically  smaller  than  Seyfert  2  ULIRGs,  which  is  rea¬ 
sonable  because  unattenuated  emission  from  unobscured  AGN s  in 
Seyfert  1  ULIRGs  can  reduce  the  EW3  3Pah  values  more  strongly 
than  attenuated  emission  from  obscured  AGNs  in  Seyfert  2 
ULIRGs.  In  Seyfert  2  ULIRGs,  the  significantly  smaller  EW3  3pah 
values  than  seen  in  starburst  galaxies  suggest  that  the  strong  dilu¬ 
tion  of  the  PAH  emission  by  the  PAH-free  continuum  produced  by 
the  AGN  is  present.  Our  L- band  spectroscopic  method  has  suc¬ 
ceeded  in  detecting  obscured  AGN  signatures  in  all  of  the  known 
obscured-AGN-possessing  Seyfert  2  ULIRGs  observed,  demon¬ 
strating  that  this  method  is  very  effective. 

Among  the  eight  Seyfert  2  ULIRGs’  nuclei,  Nagar  et  al. 
(2003)  presented  VLA  data  for  four  sources,  and  three  sources 
show  core  emission  (Table  3).  The  VLA  core  detection  rate  is 
75%  (3/4).  For  the  seven  Seyfert  1  ULIRGs’  nuclei  observed  by 
us,  five  sources  have  VLA  data  and  all  show  radio  core  emission 
(100%  =  5/5;  Table  3).  These  detection  rates  are  again  clearly 
higher  than  those  for  ULIRGs  with  no  obvious  AGN  signatures 
in  our  L-band  spectra. 

Figure  8a  plots  the  Z3 ,3pah/^ir  ratios  of  ULIRGs,  separated 
by  optical  spectral  types.  Both  non-Seyfert  and  Seyfert  ULIRGs 
are  plotted.  The  median Z3  3PAh/Tir  values  are  0.1  (LINER),  0.2 
(H  11),  0.09  (Seyfert  2),  and  <0.15  (Seyfert  1).  The  Z,3.3pah/^ir 
ratios  do  not  differ  significantly  between  non-Seyfert  and  Seyfert 
ULIRGs,  suggesting  that  both  types  of  ULIRGs  possess  de¬ 
tectable  modestly  obscured  starburst  activity  with  similar  abso¬ 
lute  luminosities. 

Figure  8  b  plots  the  rest-frame  equivalent  widths  of  the  3.3  fim 
PAH  emission  (EW3.3PAH),  separated  by  optical  spectral  types. 


In  this  plot,  Seyfert  2  ULIRGs  tend  to  show  smaller  EW3 3PAH 
values  than  non-Seyfert  ULIRGs,  while  Seyfert  1  ULIRGs  show 
the  smallest  values.  The  median  EW3 jpah  values  are  75  (LINER), 
85  (H  11),  15  (Seyfert  2),  and  <4  nm  (Seyfert  1).  Even  though  the 
absolute  luminosities  of  the  detectable  modestly  obscured  starbursts 
are  similar,  less  obscured  AGNs  can  have  a  larger  contribution  to 
the  observed L-band  flux  and  thereby  decrease  the  EW3  3PAH  values. 
The  trends  in  Figures  8a  and  8b  can  be  explained  by  a  scenario 
where  AGNs  in  Seyfert  2  ULIRGs  tend  to  be  less  obscured  by  dust 
along  our  line  of  sight  than  AGNs  in  non-Seyfert  ULIRGs.  In  fact, 
while  many  non-Seyfert  ULIRGs  show  clear  absorption  features  at 
3.1  ym  by  H20  ice  and  at  3.4  ym  by  bare  carbonaceous  dust  grains 
(Figs.  1,  2,  3),  no  such  absorption  features  are  found  in  the  L-band 
spectra  of  Seyfert  2  ULIRGs,  with  the  exception  of  IRAS  12072  - 
0444N  (Fig.  4),  supporting  the  less  obscuration  scenario  for  AGNs 
in  Seyfert  2  ULIRGs.  We  argue  that  buried  AGNs  in  non-Seyfert 
ULIRGs  and  optically  visible  AGNs  (obscured  by  torus-shaped 
dust)  in  Seyfert  2  ULIRGs  are  different  not  only  in  the  dust  ge¬ 
ometry  but  also  in  the  dust  column  density  along  our  line  of  sight,  in 
that  the  dust  columns  toward  the  AGNs  in  non-Seyfert  ULIRGs 
are  much  higher.  Since  the  dust  covering  factor  around  buried 
AGNs  in  non-Seyfert  ULIRGs  (almost  all  directions)  is  also  larger 
than  AGNs  in  Seyfert  2  ULIRGs  (torus  shaped),  the  total  amount 
of  nuclear  dust  must  be  larger  in  the  former.  The  optical  (non)- 
detectability  of  Seyfert  signatures  in  ULIRGs  depends  not  on  the 
presence  or  absence  of  powerful  AGNs,  but  largely  on  the  mass  of 
gas  and  dust  in  the  nuclei. 

5.8.  Dependence  on  Far-Infrared  Colors 

Based  on  the  IRAS  25  ^m-to-60  flux  ratio  (fislfbo), 
ULIRGs  are  divided  into  warm  (>0.2)  and  cool  (<0.2)  sources 
(Sanders  et  al.  1 988).  The  majority  of  non-Seyfert  ULIRGs  show 
cool  far-infrared  colors  (Table  1).  Optical  spectroscopy  of  gal¬ 
axies  with  warm  far-infrared  colors  often  reveals  Seyfert  sig¬ 
natures  (de  Grijp  et  al.  1987;  Keel  et  al.  2005),  while  starburst 
galaxies  usually  show  cool  far-infrared  colors.  The  warm  colors 
of  galaxies  showing  optical  Seyfert  signatures  are  usually  inter¬ 
preted  as  emission  from  dust  in  the  close  vicinity  of  an  AGN, 
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Fig.  8. — {a)  Observed  nuclear  3.3  gm  PAH-to-infrared  luminosity  ratio,  separated  by  optical  spectral  types.  ULIRGs’  nuclei  classified  optically  as  LINERs,  H  n  regions, 
unclassified,  Seyfert  2,  and  Seyfert  1  ( left  to  right)  are  separated  by  the  vertical  dashed  lines.  Filled  circles :  LINER; filled  squares :  H  n; filled  stars :  optically  unclassified;/?//^ 
triangles :  Seyfert  2;  open  circles :  Seyfert  1.  Objects’  order  is  the  same  as  that  shown  in  Table  3.  The  data  point  of  IRAS  12127— 1412  (EW3.3PAH  =  0  nm)  is  plotted  at  the 
lowest  visible  part  in  this  figure,  (b)  Rest-frame  equivalent  widths  of  the  3.3  gm  PAH  emission,  separated  by  optical  spectral  types.  Symbols  are  the  same  as  in  (a). 

whose  temperature  can  be  higher  than  that  in  starbursts  because  Figure  9b  compares  the  EW3  3PAH  and  far-infrared  colors  for 

of  the  higher  radiation  density  in  the  former.  It  is  sometimes  non-Seyfert  ULIRGs.  The  three  ULIRGs  with  almost  pure 

argued  that  cool  components  of  the  far-infrared  emission  from  buried  AGNs,  IRAS  08572+3915NW,  IRAS  12127— 1412NE, 

ULIRGs  must  be  starburst  dominated  (e.g.,  Downes  &  Solomon  and  IRAS  17044+6720,  are  distributed  on  the  warm  side,  but 

1998).  However,  the  far-infrared  color  from  AGN-heated  dust  IRAS  12127—  1412NE  has  a  cool  far-infrared  color.  Thus,  there 
is  highly  dependent  on  the  column  density  of  dust  around  an  is  at  least  one  example  of  a  pure  buried  AGN  that  produces  a  cool 

AGN,  in  that  a  larger  column  density  makes  the  far-infrared  color  far-infrared  color.  The  bulk  of  ULIRGs  with  tentative  evidence 

cooler  (/25//6O  rati0  smaller).  Our  Z-band  spectra  suggest  that  of  buried  AGNs  are  also  classified  as  cool  sources,  and  their  far- 

buried  AGNs  in  non-Seyfert  ULIRGs  are  more  highly  dust  ob-  infrared  color  distribution  overlaps  with  the  colors  of  ULIRGs 

scured  and  thus  produce  cooler  far-infrared  colors  than  AGNs  in  with  no  AGN  signatures. 

Seyfert  2  ULIRGs.  We  next  compare  the  properties  of  the  two  pure  buried 

Figure  9a  compares  the  observed  3.3  //in  PAH-to-infrared  AGNs,  IRAS  08572+3915NWand  IRAS  12127-  1412NE.  IRAS 
luminosity  ratios  with  IRAS  25  /xm-to-60  pm  flux  ratios  (i.e.,  08572+3915NW  shows  only  a  3.4  /nm  dust  absorption  feature, 

far- infrared  color)  for  non-Seyfert  ULIRGs.  There  is  no  clear  with  no  detectable  3 . 1  /nm  H20  ice  absorption.  Assuming  Galactic 

trend.  If  the  detected  3.3  pm  PAH  emission  only  probes  the  dust  properties,  Ay  >110  mag  is  estimated  toward  the  inner 

modestly  obscured  energetically  insignificant  surface  starbursts  3  pm  continuum-emitting  regions,  close  to  the  central  AGN.  On 
(§5.1)  and  buried  energy  sources  are  responsible  for  the  infrared  the  other  hand,  the  bare  dust  column  of  IRAS  12127— 1412NE 

emission  at  25  and  60  pm,  no  trend  is  expected.  (75.4  =  0.35)  is  estimated  to  be  Ay  >  50  mag.  This  source  displays 


OO  0.1  0.2  0.3  0  0.1  0.2  0.3 

IRAS  25/Ltm  to  60  jum  ratio  IRAS  25/um  to  60  /um  ratio 


Fig.  9. — (a)  Observed  nuclear  3.3  gm  PAH-to-infrared  luminosity  ratio  (ordinate)  and  IRAS  25  /zm-to-60  gm  flux  ratio  (abscissa)  for  non-Seyfert  ULIRGs’  nuclei.  Filled 
circles :  LINER  ULIRGs’  nuclei ;  filled  squares :  H  n  ULIRGs’  nuclei;  filled  stars’,  optically  unclassified  ULIRGs’  nuclei.  For  several  double-nuclei  ULIRGs,  emission  from 
individual  nuclei  is  resolved  in  our  L-band  spectra,  but  not  in  the  IRAS  data.  For  these  sources,  we  assume  that  both  nuclei  have  the  same  far-infrared  colors  as  measured  with 
IRAS.  The  data  points  of  three  LINER  ULIRGs  powered  by  almost  pure  buried  AGNs,  IRAS  08572+391 5NW,  IRAS  1 2127- 14 12NE,  and  IRAS  1 7044+6720,  are  indicated 
as  “IR08572,”  “IR12127,”  and  “IR17044,”  respectively.  The  data  point  of  IRAS  12127— 1412  (EW3.3PAH  =  0  nm)  is  plotted  at  the  lowest  visible  part  in  this  figure.  ( b )  Rest- 
frame  equivalent  widths  of  the  3.3  gm  PAH  emission  (ordinate)  and  IRAS  25  /zm— to— 60  gm  flux  ratio  (abscissa)  for  non-Seyfert  ULIRGs.  Filled  stars :  ULIRGs  with  clear 
buried  AGN  signatures;  filled  squares’.  ULIRGs  with  tentative  buried  AGN  signatures;  open  circles :  ULIRGs  with  no  obvious  buried  AGN  signatures  in  our  L-band  spectra. 
The  data  points  of  IRAS  08572+3915NW,  IRAS  12127- 1412NE,  and  IRAS  17044+6720  are  indicated  in  a  similar  way  to  (a). 
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ice  absorption  (t3  i),  and  the  column  of  ice-covered  dust  is 
estimated  to  be  A  v>6  mag.  The  total  dust  column  in  IRAS 
08572+39 15NW  may  be  higher  than  IRAS  12127- 1412NE, 
and  yet  the  former  source  shows  a  wanner  infrared  color.  This 
apparently  contradicts  the  simple  prediction  that  the  far-infrared 
color  is  cooler  when  a  larger  amount  of  dust  obscures  an  AGN, 
assuming  that  the  dust  properties  are  similar.  However,  IRAS 
08572+39 15NW  and  IRAS  12127- 1412NE  are  different  in  that 
strong  ice  absorption  is  found  only  in  the  latter  source.  Ice- 
covered  dust  grains  can  make  a  cooler  far-infrared  color  (smaller 
fzs  If 60  ratio)  than  do  bare  dust  grains  because  of  the  change  of 
optical  constants  in  the  infrared  ( Dudley  et  al.  2003).  The  cooler 
color  of  IRAS  12127-1412NE  than  IRAS  08572+3915NW 
may  be  due  to  the  presence  of  a  significant  amount  of  ice. 

6.  SUMMARY 

We  presented  the  results  of  infrared  Z-band  (3-4  fi m )  slit 
spectroscopy  of  the  nuclei  of  a  large  sample  of  nearby  ULIRGs 
in  the  IRAS  1  Jy  sample.  ULIRGs  with  no  obvious  Seyfert  sig¬ 
natures  in  the  optical  (i.e.,  LINER,  H  n,  and  unclassified)  were 
mainly  observed  and  spectra  of  27  LINER,  13  H  n,  and  two 
optically  unclassified  ULIRGs’  nuclei  were  obtained.  Using  the 
3.3  //m  PAH  emission,  3.1  //m  H20  ice  absorption,  and  3.4  ft m 
carbonaceous  dust  absorption  features,  we  investigated  whether 
the  cores  of  these  ULIRGs  can  be  explained  by  very  compact 
starbursts  only  or  show  signatures  of  buried  AGNs.  Since  our 
spectra  covered  both  the  shorter  and  longer  wavelength  sides  of 
these  features,  we  were  able  to  distinguish  between  absorption- 
dominated  and  emission-dominated  sources  with  small  ambigui¬ 
ties.  These  spectra  were  compared  with  known  AGN-possessing 
ULIRGs’  nuclei  (eight  Seyfert  2  and  seven  Seyfert  1),  and  we 
found  the  following  main  conclusions: 

1.  The  3.3  fim  PAH  emission,  the  probe  of  starbursts,  was 
detected  in  all  but  two  non-Seyfert  ULIRGs.  It  is  likely  that 
modestly  obscured  (Ay  <  15  mag)  starbursts  at  the  surface  of 
ULIRGs’  cores  are  responsible  for  the  detected  PAH  emission. 
The  observed  3.3  fim  PAH-to-infrared  luminosity  ratios  were 
smaller  by  a  factor  of  2  to  >10  than  those  found  in  less  infrared- 
luminous  starbursts,  suggesting  that  the  detected  modestly  ob¬ 
scured  starbursts  are  energetically  insignificant. 

2.  Three  LINER  ULIRGs,  IRAS  08572+39 15NW,  IRAS 
12127— 1412NE,  and  IRAS  17044+6720,  showed  3.3  fim  PAH 
emission  equivalent  widths  substantially  smaller  than  typical 
values  found  in  starburst  galaxies.  These  ULIRGs  were  classi¬ 
fied  as  sources  dominated  by  almost  pure  buried  AGNs,  with  no 
detectable  or  very  weak  starburst  activity. 

3 .  Besides  these  three  ULIRGs,  the  absolute  optical  depths  of 
the  3.1  /jm  FLO  ice  absorption  feature  were  used  to  distinguish 
whether  the  energy  sources  at  the  cores  of  ULIRGs  are  more 
centrally  concentrated  than  the  surrounding  dust  (as  is  expected 
for  a  buried  AGN)  or  are  spatially  well  mixed  with  dust  (a  nor¬ 
mal  starburst).  The  large  optical  depths  in  14  LINER,  three  H  n, 
and  two  unclassified  ULIRGs’  nuclei  suggested  that  the  energy 
sources  are  more  centrally  concentrated  than  the  surrounding  dust, 
unless  a  large  amount  of  foreground  ice-covered  dust  is  present 
in  the  interstellar  medium  of  the  host  galaxies,  which  we  consider 
unlikely.  These  ULIRGs  were  classified  as  showing  tentative  evi¬ 
dence  for  powerful  buried  AGNs. 

4.  Based  on  the  optical  depths  of  the  3.4  fim  carbonaceous 
dust  absorption  feature,  centrally  concentrated  energy  sources 
were  suggested  in  the  two  LINER  ULIRGs,  IRAS  08572+ 
3915NW  and  IRAS  12127-  1412NE.  Both  of  them  had  already 
been  classified  as  buried  AGNs  from  the  3.3  fim  PAH  equivalent 


widths,  and  thus  powerful  buried  AGNs  were  consistently  iden¬ 
tified  with  two  independent  methods. 

5.  In  total,  17  out  of  27  (63%)  LINER  ULIRGs’  nuclei  and  3 
out  of  13  H  ii  ULIRGs’  nuclei  (23%)  showed  some  evidence  for 
buried  AGNs.  The  higher  detection  rate  of  buried  AGNs  in 
LINER  ULIRGs  than  in  H  n  ULIRGs  was  similar  to  the  higher 
detection  rate  of  compact  radio  cores  (also  a  signature  of  AGNs) 
in  the  former. 

6.  Z-band  spectra  of  eight  Seyfert  2  and  seven  Seyfert  1 
ULIRGs’  nuclei  were  obtained,  and  the  3 .3  fim  PAH  emission  was 
detected  in  seven  (out  of  eight)  and  five  (out  of  seven)  sources, 
respectively,  suggesting  that  these  Seyfert  ULIRGs  also  contain 
detectable  starburst  activity. 

7.  Based  on  the  equivalent  widths  of  the  3.3  fim  PAH  emis¬ 
sion,  our  Z-band  spectra  successfully  detected  AGN  signatures 
in  all  of  the  observed  Seyfert  2  ULIRGs  (known  obscured- AGN- 
possessing  ULIRGs),  demonstrating  that  our  /.-band  spectroscopic 
method  is  indeed  effective  in  revealing  signatures  of  obscured 
AGNs. 

8.  The  detection  rate  of  VLA  radio  cores  was  higher  in 
ULIRGs  that  display  evidence  of  buried  AGNs  in  the  L  band. 
This  suggested  that  our  Z-band  spectroscopic  method  is  reliable, 
if  not  perfect. 

9.  Among  eight  observed  Seyfert  2  ULIRGs,  only  IRAS 
12072— 0444N  showed  a  clear  3.4  fim  carbonaceous  dust  ab¬ 
sorption  feature.  The  remaining  Seyfert  2  ULIRGs  showed 
neither  detectable  3.4  fim  nor  3.1  fim  absorption  features. 

10.  There  was  no  systematic  difference  in  the  3.3  fim  PAH- 
to-infrared  luminosity  ratios  between  Seyfert  and  non-Seyfert 
ULIRGs.  However,  Seyfert  1  and  Seyfert  2  ULIRGs  tended  to 
show  smaller  rest-frame  equivalent  widths  of  the  3.3  fim  PAH 
emission  than  non-Seyfert  ULIRGs,  with  Seyfert  1  ULIRGs 
(unobscured  AGNs)  exhibiting  the  smallest  values.  The  optical 
depths  of  absorption  features  in  the  Z  band  were  also  generally 
smaller  in  Seyfert  2  ULIRGs  than  in  non-Seyfert  ULIRGs.  These 
results  were  naturally  explained  if  AGN  s  in  non-Seyfert  ULIRGs 
are  more  highly  obscured  by  dust  along  our  sight  line  than  the 
AGNs  in  Seyfert  2  ULIRGs. 

11.  Evidence  for  buried  AGNs  was  found  in  non-Seyfert 
ULIRGs  with  both  warm  and  cool  far-infrared  colors.  The 
LINER  ULIRG  dominated  by  a  pure  buried  AGN,  IRAS  12127— 
1412NE,  has  a  cool  far- infrared  color.  The  cool  far-infrared  color 
in  this  pure  buried  AGN  could  be  the  consequence  of  a  large 
dust  column  density  around  an  AGN,  including  ice-covered 
dust,  as  is  implied  from  strong  dust  absorption  features  at  3.1 
and  3.4  fi,m. 

12.  Our  Z-band  spectroscopy  revealed  the  presence  of  AGNs 
in  many  non-Seyfert  ULIRGs,  and  these  AGNs  are  thought  to  be 
obscured  (at  least  opaque  to  the  AGN’s  ionizing  radiation)  along 
virtually  all  sight  lines  so  as  to  be  optically  elusive  (Maiolino 
et  al.  2003).  This  dust  geometry  is  fundamentally  different  from 
the  AGNs  in  Seyfert  2  ULIRGs,  which  are  obscured  by  dust  in 
a  torus  geometry.  Together  with  the  above  finding  that  dust 
column  along  our  line  of  sight  is  higher  in  AGNs  in  non-Seyfert 
ULIRGs,  it  was  implied  that  a  larger  amount  of  dust  is  con¬ 
centrated  surrounding  the  AGNs  at  the  nuclei  of  non-Seyfert 
ULIRGs  than  Seyfert  ULIRGs.  We  argued  that  the  (non)- 
detection  of  Seyfert  signatures  in  the  optical  is  highly  depen¬ 
dent  on  the  total  amount  of  dust  and  its  covering  factor  around  a 
central  AGN  and  that  many  non-Seyfert  ULIRGs  may  contain 
powerful  buried  AGNs.  Other  types  of  careful  investigations  at 
the  wavelengths  of  lower  dust  extinction  are  important  to  find 
new  buried  AGNs  that  may  even  be  undetectable  in  our  Z-band 
spectroscopic  method. 
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